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Effects  of  temperature  on  the  developmental  rates  of  Anticarsia 
gemmatalis  (Hubner)  and  an  associated  nucleopolyhedrosis  virus  were 
investigated  through  a  series  of  experiments  at  various  constant 
temperatures  on  artificial  diet.   Photoperiod  in  the  A.  gemmatalis 
experiments  was  12:12  photophase:scotophase;  photoperiod  in  the 
nucleopolyhedrosis  virus  experiment  was  14:10  photophase:scotophase. 
Humidity  was  not  controlled;  however,  relative  humidity  was  never 
below  50%.   Parameters  for  a  general  poikilo thermic  developmental 
model  were  estimated  for  both  organisms.   Validation  data  were 
collected  from  "laboratory  reared"  and  naturally  occurring  A. 
gemmatalis  larvae  held  at  constant  and  variable  temperatures  in  both 
laboratory  and  field.   Validation  data  for  the  nucleopolyhedrosis 
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virus  system  were  collected  from  "laboratory  reared"  larvae  held 
under  controlled  variable  temperatures. 

Survival  of  A^  gemmatalis  cohorts  held  at  various  constant 
temperatures  changed  little  between  21.1  and  32.2  C.   Maximum 
survival  occurred  at  23.9°C.  Mortality  increased  substantially  at 
temperature  extremes  (15.6,  18.3,  and  37.8°C)  with  maximum  mortality 
occurring  at  37.8°C.   In  general,  mortality  occurred  at  the 
metamorphosis  into  the  pupal  stage  and  was  accompanied  with  definitive 
signs  of  stress  induced  septicemia.   Developmental  times  consistently 
decreased  with  increasing  temperature  up  to  32.2  C.   Higher  tempera- 
tures were  associated  with  increasing  developmental  times. 

Proportion  of  the  total  developmental  time  spent  in  any  one  size 
class  ranged  from  .32  to  .42  for  oviposition  to  1/2"  larvae,  .12  to 
.17  for  1/2"  larvae  to  1"  larvae,  .11  to  .17  for  1"  larvae  to 
pupation,  and  .27  to  .38  for  pupation  to  adult  eclosion.   Statistical 
analysis  indicated  that  these  proportions  did  not  remain  constant  over 
all  temperatures  tested.  The  ratio  of  standard  deviation  to  mean 
developmental  rate  was  computed  for  all  constant  temperatures.   These 
ratios  did  not  appear  to  remain  constant.   Comparisons  of  develop- 
mental times  of  "laboratory  reared"  and  insects  taken  from  naturally 
occurring  populations  were  made  under  three  temperature  regimes . 
Statistical  analysis  demonstrated  significant  difference  in  6  of  12 
trials.   Maximum  difference  was  2.1  days,  while  minimum  difference 
was  0  days. 

Emergence  patterns  predicted  by  the  developmental  model  were 

compared  with  observed  emergence  of  A^  gemmatalis  in  laboratory  and 

field.   In  general,  these  patterns  closely  agreed.   Reasons  for 

departure  from  some  model  predictions  were  discussed. 
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Developmental  times  of  the  nucleopolyhedrosis  virus  of  A. 
gemmatalis  decreased  with  increasing  temperature  to  30.0  C.   Higher 
temperatures  were  associated  with  increasing  developmental  times. 
Viral  replication  was  inhibited  completely  at  40  C.   Temperature 
and  dose  level  were  shown  to  have  significant  influence  on  viral 
developmental  times  and  upon  developmental  times  of  surviving  A^ 
gemmatalis  larvae.   Proportion  of  the  total  developmental  time  of 
A.  gemmatalis  spent  in  the  viral  developmental  stage  was  shown  by 
statistical  analysis  not  to  remain  constant  over  various  temperatures. 

Occurrence  of  observed  mortality  and  mortality  predicted  by  the 
model  were  compared  at  four  variable  temperature  regimes.   In  general, 
observed  and  predicted  time  of  mortality  agreed.   Reasons  for 
departure  from  some  model  predictions  were  discussed. 
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CHAPTER  I 
INTRODUCTION 


Soybeans  were  first  grown  in  Florida  on  an  experimental  basis  in 
ca.  1925.   By  the  late  I9A0's,  soybean  became  economically  important 
in  Florida  (Hinson  1967) .   Production  of  soybean  has  increased  in 
Florida  and  throughout  the  United  States;  subsequently,  1965  acreage 
of  soybean  in  Florida  was  ca.  30,400  ha.  (76,000  acres). 

The  Florida  soybean  agroecosystem  presently  (1980)  consists  of 
ca.  113,600  ha.  (334,000  acres),  grown  primarily  in  northcentral  and 
northwest  Florida  (Koehn  19  78).   Soybeans  are  attacked  by  a  variety 
of  foliage  and  pod  feeding  insects.   The  Florida  Insect  Control  Guide 
(Florida  Agricultural  Experiment  Station)  lists  recommended  controls 
for  selected  soybean  pod  feeders:   corn  earworm,  Heliothis  zea 
(Boddie) ;  southern  green  stinkbug,  Nezara  viridula  (L.);  and  foliage 
feeders:   velvetbean  caterpillar,  Anticarsia  gemmatalis  (Hubner) ; 
mexican  bean  beetle,  Epilachna  varivestis  Mulsant;  cabbage  looper, 
Trichoplusia  ni  (Hubner);  soybean  looper,  Pseudoplusia  includens 
(Walker) .   The  most  damaging  foliage  feeder  typically  encountered  in 
Florida  is  the  velvetbean  caterpillar  (VBC)  (Fla.  Agric.  Exp.  Sta. 
1974,  Strayer  and  Greene  1974). 

Nickles  (1926)  was  first  to  report  infestation  of  soybean  by  VBC 
in  Florida,  and  Watson  (1916)  reported  the  presence  of  VBC  in  Florida 
as  a  major  pest  of  velvetbeans,  Stizolobium  deeringianum  Bort.   As 
acreage  of  soybean  production  increased,  so  increased  problems 


caused  by  VBC.   Strayer  (1973)  estimated  that  the  annual  cost  of 
controlling  VBC  could  reach  $1.2  million. 

At  present,  the  VBC  is  not  considered  to  have  the  ability  to 
overwinter  in  the  soybean  production  region  of  Florida  (Buschman 
et  al.  1977,  Stimac  unpublished  data).   Watson  (1916)  indicated  that 
the  VBC  overwintered  in  south  Florida  and  flew  northward  each  year. 
However,  VBC  populations  may  be  able  to  overwinter  as  pupae  much 
farther  north,  at  least  in  mild  years  (Buschman  et  al.  1977).   Stimac 
(unpublished  data)  found  that  too  few  pupae  could  survive  in  north 
Florida  to  produce  economic  infestations  the  next  season. 

VBC  populations  are  subject  to  epizootics  of  the  fungal  pathogen 
Nomuraea  rileyi  (Farlow)  Sampson  (Allen  et  al.  1971).   This  entomo- 
pathogen  has  been  observed  for  many  years  to  be  a  significant  agent  of 
VBC  mortality  during  late  season  (Watson  1916,  Kuitert  1967).   A 
nucleopolyhedrosis  virus  (VBC-NPV)  imported  from  Brazil  (Allen  and 
Knell  1977)  has  been  demonstrated  to  have  potential  use  in  a 
management  program  directed  at  the  VBC.   However,  precise  definition 
of  the  VBC  life  system  must  precede  attempts  at  VBC  management. 

Rabb  (1970,  p.  2)  defines  a  pest  management  system  as 
"...  the  reduction  of  pest  problems  by  actions  selected  after  the 
life  systems  of  the  pest  are  understood  and  the  ecologic  as  well  as 
the  economic  consequences  of  these  actions  have  been  predicted,  as 
accurately  as  possible  to  be  in  the  best  interest  of  mankind." 
Theoretically,  pest  management  is  not  composed  of  a  particular  tactic 
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but  a  combination  of  compatible  tactics  decided  upon  after  relevant 
aspects  of  the  pest  life  system  are  understood.   To  insure  optimum 
results,  the  compatibility  of  combinations  of  individual  tactics,  in 
the  context  of  the  target  system,  must  be  understood. 

Berryman  and  Pienaar  (1974)  state  that  the  acceptance  of  the 
need  for  complex  crop-pest  information,  such  as  that  suggested  by  Rabb 
(1970)  necessitates  the  formulation  of  a  strategic  approach  and  a 
formalization  of  decision-making  policy  for  pest  management.   However, 
the  acceptance  of  these  concepts  leads  to  the  realization  that  the 
outcomes  of  events  in  highly  complex  natural  systems  are  difficult 
to  predict.   These  events  are  often  counterintuitive;  i.e.,  the 
actions  of  the  system  do  not  lead  to  what  one  would  expect  (Watt 
1970).   The  fundamental  concept  of  systems  theory  is  hierarchic 
order.   Stimac  and  Overton  (1979)  suggest  the  use  of  systems  science 
as  a  method  of  deciphering  complex  systems.   They  state  that  while 
systems  science  will  provide  the  techniques  for  construction  and 
analysis  of  pest-crop  systems,  discipline  scientist  must  provide 
information  necessary  to  identify  and  quantify  relevant  parameters. 

Figures  1-3  (Stimac  and  Barfield  1980)  detail  a  conceptual  model 
of  the  Florida  soybean  agronomic  system.   The  hierarchical  nature  of 
this  system  provides  for  changes  in  one  subsystem  to  affect  other 
subsystems,  and  focuses  on  the  importance  of  understanding  relevant 
system  processes  (e.g.,  development,  mortality,  and  movement).   Typical 
life  history  of  any  single  insect  pest  in  the  system  is  depicted  in 
Figure  3.   Emphasis  at  the  pest  population  level  is  on  processes  of 
development,  mortality,  consumption,  and  oviposition. 
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Estimates  of  VBC  soybean  foliage  consumption  and  ovipositional 
rates  have  been  provided  by  Moscardi  (1979).   Luna  (1979)  estimated 
VBC  larval  stage  mortality.  However,  information  must  exist  on  the 
temperature-dependent  development  periods  for  at  least  relevant 
"size  classes"  of  pests  prior  to  any  management  strategies  being 
invoked.   The  lack  of  such  information  on  the  key  insect  pest  (VBC) 
of  Florida  soybean  has  stimulated  present  investigations  toward  the 
construction  of  a  VBC  developmental  model  specifically  designed  to 
provide  estimates  of  developmental  times  as  a  function  of  temperature. 

The  simultaneous  construction  of  a  VBC-NPV  developmental  model 
was  necessary  to  estimate  NPV- induced  insect  mortality  as  a  function 
of  temperature.   Both  models  are  needed  to  provide  information 
critical  to  the  management  of  VBC  within  the  soybean  agroecosystem. 


CHAPTER  II 
LITERATURE  REVIEW 


Velvetbean  Caterpillar 

A  vast  body  of  literature  has  been  compiled  relative  to  the 
velvetbean  caterpillar  (VBC) ,  Anticarsia  gemmatalis  (Hubner) .   Life 
history  data  under  natural  (Watson  1915;  1916,  Douglas  1930,  Hinds 
1930,  Hinds  and  Osterberger  1931,  Ellisor  1942)  and  controlled 
conditions  (Leppla  et  al.  1977)  have  been  documented.   Mating  and 
ovipositional  behavior  were  reported  by  Greene  et  al.  (1973),  and 
behavior  associated  with  mating  and  pheromone  production  by  Johnson 
(1976).   Reid  (1975)  and  Nickle  (1977)  reported  VBC  consumption  of 
soybean  and  pf.inut  leaves,  respectively.   Strayer  (1973)  developed  a 
sequential  sampling  plan  for  VBC  larvae  on  soybean.   Comprehensive 
lists  of  VBC  host  plants,  parasites,  predators  and  pathogens  can  be 
found  in  Moscardi  (1979). 

Literature  concerning  the  effect  of  temperature  on  the  VBC  is 
relatively  scarce  and  disjimct.  Nickle  (1977)  reared  VBC  on  peanut, 
Arachis  hypogaea  L.  cv.  "Florunner,"  leaves  at  constant  temperatures 
of  21.2,  23.8,  26.7,  29.4,  32,2,  and  35  +  1°C.   Larvae  required 
averages  of  29.7,  25.7,  24.2,  22.5,  21.2,  and  20.0  days,  respectively, 
to  develop  from  eggs  to  pupae.   Mean  pupal  duration  at  these  same 
constant  temperatures  required  14.3,  12.5,  10.0,  8.0,  8.4,  and  6.7 
days,  respectively.   Reid  (1975)  reared  insects  on  soybean  foliage 
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at  constant  temperatures  of  15.6,  18.3,  23.9,  and  29.4  +  1.5°C. 
Duration  of  the  pupal  stage  was  reported  at  19.4  and  7.7  days  for 
temperatures  18.3  and  29.4  +  1,5  C,  respectively. 

Leppla  et  al.  (1977)  determined  the  larval  survival,  time  to 
pupation,  time  to  adult  eclosion,  adult  survival  and  oviposition  of 
VBC  at  21.1,  26.7,  and  32.2  +  1  C.   Mean  larval  development  times 
(days)  for  males : females  were  30.76:39.8,  15.5:15.74,  and  12.5:12.4. 

Moscardi  (1979)  reported  the  effect  of  temperature  on  egg  lay, 
egg  hatch  and  longevity  for  VBC  adults.   Using  constant  temperatures 
of  21.1,  23.9,  26.7,  29.4,  and  32.2  +  1  C,  the  mean  eggs/female  were 
482,  732,  842,  and  713,  respectively.   Mean  percent  egg  hatch  was 
68.6,  77.4,  82.3,  78.5,  and  80.8,  respectively,  for  these  same 
temperatures.   Mean  longevity  for  adult  females  (days)  was  24.8,  21.9, 
18.0,  15.5,  and  11.2,  respectively. 

Poikilothermy 

The  effect  of  temperature  on  living  organisms  has  been  the 
subject  of  many  studies;  thus  a  bibliography  of  the  subject  would  be 
tedious.   An  extremely  thorough  review  of  this  subject  can  be  found 
(Precht  et  al.  1973)  which  covers  plants,  microbes,  invertebrate  and 
vertebrate  animals.   Concepts  from  the  subcellular  to  the  population 
levels  are  explained  in  detail,  and  most  of  the  existing  literature 
dealing  with  these  areas  is  utilized.   Consequently,  only  an  overview 
of  how  temperature  affects  and  interacts  with  insects  will  be 
presented  here. 

Poikilothermic  temperature  regulation  in  organisms  is  described 
by  Jankowsky  (1973)  as  follows:   the  body  temperature  of  poikilotherms 
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changes  in  the  same  direction  as  ambient  temperature  but  is  rarely 
the  same.   Most  organisms  on  earth  are  poikilo therms:   plants, 
invertebrates,  fish,  amphibians,  and  reptiles.   Birds  (Aves)  and 
mammals  (Mammalia)  are  true  homeotherms.   However,  the  distinction 
between  the  types  of  temperature  regulation  has  not  always  been  clear. 
Poikilothermy  can  be  divided  into  two  subtypes,  heterothermy  and 
ectothermy.   The  former  refers  to. the  increase  of  body  temperature  by 
endogenous  heat  production;  the  latter  refers  to  modification  in 
external  heat  input  through  behavioral  traits  of  the  organism 
(Jankowsky  1973) . 

Invertebrates,  in  general,  have  a  very  low  body  weight  and  thus 
a  small  capacity  to  produce  and  retain  endogenous  heat .   Thus ,  some 
organisms  (e.g.,  insects)  have  behavioral  characteristics  conducive 
to  heat  generation.   Raising  the  basal  metabolism  by  shivering  or 
flapping  the  wings  can  be  used  to  warm  insect  flight  muscles 
(Hannegan  and  Heath  1970).   Many  insects  may  require  solar  heating 
before  becoming  active  (Digby  1955) .   Raising  the  body  temperature 
through  external  radiation  may  be  accomplished  through  such  means 
as  coloration  (arctic  Lepidoptera  are  often  dark  in  color)  or  by 
increasing  body  size,  as  in  Orthoptera,  Hymenoptera  and  Diptera 
(Digby  1955).   Changes  in  body  surface,  rotation  of  the  body  axis, 
spreading  of  wings  and  movement  to  and  from  warm  locations  may  also 
serve  to  control  body  temperature  (Flitters  1968) .   Social  insects 
can  control  the  temperature  of  the  colony  by  insulating  the  nest,  as 
in  some  ants;  production  of  heat  through  the  pumping  of  the  abdomen, 
as  in  wasps  and  hornets;  basking  in  the  sun;  or  movement  to  the 
interior  of  the  nest  (Jankowsky  1973) . 
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Temperature  Effects  on  Insects 

The  effect  of  temperature  on  insects  is  so  pronounced  that  most 
scientific  papers  dealing  with  some  phase  of  insect  life  history  will 
have  at  least  cursory  mention  of  temperature  effects.   Wide-spread 
recognition  of  the  importance  of  temperature  bodes  both  good  and  bad  - 
for  the  entomological  sciences.   A  great  majority  of  temperature- 
related  insect  information  has  been  gathered  only  within  narrow 
temperature  ranges  (Hight  et  al.  1972,  Leppla  et  al.  1977,  LeCato 
and  Pienkowski  1972) . 

Temperature  effects  on  insects  can  be  direct  or  indirect, 
physiological  or  behavioral,  straight-forward  or  hidden.   The 
following  paragraph  depicts  some  ways  by  which  temperature  can 
affect  insects.   Again,  for  a  review  of  mechanistic  details,  the 
reader  is  referred  to  Precht  et  al.  (1973). 

Temperature  may  affect  a  variety  of  insect  behaviors.   Allen  and 
Gonzalez  (1975),  working  with  the  parasite,  Trichogramma  pretiosum 
Riley,  and  its  host,  the  cabbage  looper,  found  that  the  functional 
response  of  the  parasite  to  dynamic  host  density  changes  as  a  function 
of  temperature.   This  behavior  was  mediated  by  the  fact  that  low 
temperature  provided  a  "refuge"  for  the  host  due  to  reduced  search  by 
the  parasite.   High  temperatures  did  not  provide  the  "refuge."  Fewer 
hosts  being  attacked  at  higher  temperatures  were  attributed  to  the 
unavailability  of  hosts  due  to  thermal  mortality.   Biever  (1972) 
found  that  search  behavior  of  the  parasites  Trichogramma  minutum 
Riley,  T^   evanescens  Westwood,  T^  semifumatum  (Perkins) ,  and  T^  £p, 
varied  with  temperature.   Increases  in  experimental  temperatures 
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(ET  =  20,  25,  30,  35,  and  40  +  1  c)  resulted  in  increased  searching 
time;  at  40  C  search  rate  declined.   Crocker  et  al ,  (1975)  studied 
the  effect  of  temperature  on  predator  behavior  and  found  that  a 
greater  number  of  prey/day  were  consumed  as  temperature  increased 
(ET  =  20,  25,  30,  and  35+1  C) .   Here,  number  of  prey  consumed 
declined  at  35  C.   This  relationship  lead  to  the  suspicion  that  a 
sample  of  predators  in  any  particular  environment  did  not  necessarily 
provide  a  suitable  index  as  to  the  effect  of  predation  on  a  host 
population,   DeRozari  et  al.  (1977)  discovered  that  temperature 
indirectly  affected  the  mating  behavior  of  the  European  corn  borer, 
Ostrinia  nubilalis  (Hubner) ,  which  would  not  mate  unless  liquid  water 
was  present.   Temperature  has  an  effect  on  this  behavior  because  it 
is  temperature  which  mediates  dew  formation. 

Distribution  of  insect  populations  are  affected  by  temperature. 
Areas  available  to  populations  of  Drosophila  inornata  Malloch  in 
Australia  are  temperature  limited.   Temperatures  12°C  or  below,  and 
high  temperature/humidity  combinations  of  20-30°C/90%  RH  exclude 
successful  exploitation  by  these  insects  (Parsons  1975). 

Insect  morphology  and  physiology  are  affected  by  temperature. 
Zenner-Polania  and  Helgesen  (1973),  working  with  the  omnivorous  leaf 
roller,  Platynota  stultana  (Washingham) ,  found  that  within  a  cohort, 
the  proportion  of  insects  developing  through  6  instars  decreased  as 
temperature  increased  (ET  =  15,  20,  25,  30,  and  35  +  1°C) .   Mean  head 
capsule  width  decreased  at  temperatures  above  25°C.   Tuskes  and  Atkins 
(1973)  produced  summer  or  winter  coloration  respectively  in  adults 
of  the  alfalfa  caterpillar,  Colias  eurytheme  Boisduval,  by  subjecting 
juveniles  of  the  previous  generation  to  either  high  (26.6  or  32.2 
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+  .47  C)  or  low  temperatures  (18.3  +  .47  C) .   Intermediate  temperatures 
(21.2+  .47  C)  produced  mixed  proportions  of  summer  and  winter  colored 
adults.   Gorsuch  and  Karandinos  (1974)  found  that  the  daily  emergence 
of  the  lesser  peachtree  borer,  Synanthedon  pictipes  (Grote  and  Robin- 
son) ,  was  related  to  the  rate  of  temperature  increase  rather  than  a 
temperature  threshold.   Kalpage  and  Brust  (1974)  reported  that  high 
temperatures  can  terminate  the  embryonic  diapause  of  Aedes  atropalpus 
(Coquillett) .   Thurston  (1972)  found  that  induction  of  diapause  in 
the  tobacco  hornworm,  Manduca  sexta  (L.),  could  be  mediated  by  either 
a  shorter  photoperiod  or  a  lower  temperature. 

Insect  survival  is  affected  by  temperature.   Temperature  was 
found  to  moderate  the  effect  of  irradiation  on  survival  of  lesser 
grain  borer,  Rhyzopertha  dominica  (F.).   Elevation  of  temperature 
(ET  =  19,  25,  30,  and  35+1  C)  reduced  the  percent  and  duration  of 
survival  after  irradiation,  even  though  control  individuals  lived 
longer  and  had  higher  survival  at  30  C  than  at  all  other  temperatures 
(Singh  and  Liles  1972) .   The  effectiveness  of  several  insecticides 
has  been  shown  to  change  with  temperature.   Meisch  et  al.  (1972) 
tested  the  toxicity  of  azinphosmethyl  and  malathion  on  the  boll 
weevil,  Anthonomus  grandis  Boheman,  at  various  temperatures. 
Azinphosmethyl  was  found  to  be  more  toxic  at  warmer  temperatures 
(26.7  C) ,  while  malathion  was  more  toxic  at  cooler  temperatures 
(16.7  C) .   lordanou  and  Watters  (1969)  found  that  malathion  and 
bromophos  had  positive  temperature  coefficients  while  DDT  had  a 
negative  temperature  coefficient.   The  effect  of  high  temperature  on 
survival  of  the  sugarcane  borer,  Diatraea  saccharalis  (F.),  was 
investigated  by  Miskimen  (1973).   Constant  high  temperatures  were 
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found  to  be  more  detrimental  than  fluctuating  temperatures.   Insects    ,* 
can  be  sterilized  by  high  temperatures.   Cheng  (1972)  found  that  5   '  '' 
days  at  35  C,  during  the  pupal  stage,  sterilized  males  of  Eucosma    .  '  k.; 


schistaceana  Snellen.   Sterilization  was  expressed  by  the  inability 
of  males  to  transfer  sperm  and  reduced  male  competitiveness. 

Low  temperature  survival  of  the  Gypsy  moth,  Porthetria  dispar 
(L.),  was  shown  to  be  mediated  by  snowfall  (Leonard  1972).   Height  of 
egg  mass  above  ground  level  was  critical  for  production  of  the  next 
generation  because  snow  provided  insulation  from  cold.   Eggs  pro- 
tected by  snow  survived  to  -32.2°C.   David  et  al.  (1977)  demonstrated 
that  acclimation  exposure  of  gypsy  moth  larvae  to  cool  temperatures, 
before  exposure  to  extreme  cold,  increased  survival. 

Host  plant  seclection  is  another  insect  characteristic  affected 
by  temperature.   Wood  and  Starks  (1972)  found  that  greenbugs, 
Schizaphis  graminum  (Rondiani) ,  live  longer  on  preferred  hosts. 
They  also  found  that  greenbug  reproduction  is  maximum  at  different 
temperatures,  dependent  upon  host  plant  nutrition.   Increased  temper- 
ature reduced  fecundity,  length  of  lifespan,  and  reproductive  periods 
Starks  et  al  (1973)  found  that  temperature  affected  preference  for 
host  plant  (varieties  of  Sorghum  bicolor  (L.)  Moench)  by  biotypes  of 
S.  graminum.   Biotype  B  showed  an  increased  preference  for  "Deer"  and 
"Piper"  varieties  and  decreased  preference  for  "OK-8"  as  temperature 
increased.   However,  preferences  of  Biotype  C  were  not  affected  by 
temperature.   Sosa  and  Foster  (1976)  demonstrated  that  infestations 
of  the  Hessian  fly,  Mayetiola  destructor  (Say),  increased  with 
temperature  in  four  cultivars  tested. 
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Insect  reproductive  biology  is  influenced  by  temperature.   LeCato 
and  Pienkowski  (1972)  found  that  the  reproductive  efficiency  of  the 
alfalfa  weevil,  Hypera  postica  (Gyllenhal) ,  was  greatest  at  inter- 
mediate temperatures  (23  C)  while  higher  temperature  (30  C)  was 
unfavorable  for  oviposition.   Insects  held  under  a  square  wave 
temperature  regime  (12h:12h,  4.5:15.5  C)  were  more  productive  than 
insects  at  all  constant  temperatures  (10,  23,  and  30  C) .   Hennberry 
et  al.  (1977)  found  that  the  pink  bollworm,  Pectinophora  gossypiella 
(Saunders) ,  had  an  optimum  temperature  range  for  reproduction 
(27-32°C).   Temperatures  above  this  range  reduced  longevity, 
fecundity,  and  egg  variability.   Males  reared  at  35  C  mated  as 
frequently  as  other  males,  but  in  85%  of  the  matings,  no  sperm  were 
transferred.   Dunbar  and  Bacon  (1972)  found  that  oviposition  by 
Geocoris  Spp.  could  be  halted  by  either  low  (23.9  C)  or  high  (35.0  C) 
temperatures.   Bryan  et  al.  (1969)  demonstrated  that  production  of 
progeny  by  Lespesia  archippivora  (Riley)  was  optimum  at  30  C  and 
decreased  above  or  below  that  temperature.   Hagley  (1972)  found  that 
high  temperature  compounded  the  effect  of  female  age  on  percent  egg 
hatch.   Amman  (19  72)  found  that  oviposition  by  female  Dendroctonus 
ponderosae  Hopkins  increased  with  temperature  (ET  =  10,  15,  and  20  C) . 
Production  of  progeny  and  longevity  of  Eucelatoria  sp  had  temperature 
optima  above  which  efficiency  decreased  (Bryan  et  al.  1972).   The 
production  of  sexuals  in  Therioaphis  maculata  (Buckton)  and  T^  riehmi 
(Borner)  was  shown  to  be  influenced  by  temperature  (Schalk  and 
Manglitz  1972). 
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Developmental  Process 

Animal  development  takes  place  within  a  relatively  narrow  range 
of  temperatures.   Description  of  this  range  for  a  particular 
organism  involves  measurement  of  the  time  required  to  change  from 
one  life  stage  to  another  at  a  series  of  constant  temperatures.   In 
such  cases,  one  is  dealing  with  a  fully  adapted  organism  and,  as 
such,  any  adaptation  phenomena  present  are  included  within  the 
measurements  (Laudien  1973).  .  ,-   . 

The  relation  of  temperature  to  development  can  be  expressed  by 
a  general  temperature-development  profile  (Figure  4) .   The  lower 
portion  of  the  profile  is  nonlinear  and  contains  the  "critical  point" 
(developmental  zero);  the  middle  portion  is  linear;  and  the  upper 
portion  is  nonlinear,  where  rates  again  begin  to  decrease.   The 
lower  end  of  a  poikilo therm's  temperature  range  is  characterized  by 
survival  and  extremely  slow  development  (Laudien  1973).   Determining 
this  temperature  threshold  experimentally  is  difficult.   Estimating 
the  "critical  point"  with  the  temperature  at  which  development  can 
first  be  measured  was  suggested  by  Shelford  (1927).   However,  this 
approach  becomes  directly  dependent  upon  the  precision  of  the 
observer.   Pearis  (1914,  p.  11)  stated  that  this  "critical  point"  was 
"at  or  near  the  point  where  the  reciprocal  curve  for  the  time  factor 
intersects  the  temperature  axis."  However,  Pearis  suggested  a  linear 
temperature-development  rate  relationship.   Developmental  rate 
(1/time)  in  this  region  of  the  curve  does  not  decrease  in  a  linear 
fashion  but  approaches  zero  asymptotically,  as  with  the  corn  earworm 
(Mangnat  and  Apple  1966).   This  critical  point  will  vary  among  species, 
races,  or  even  stages  of  development  (Sanderson  1910). 
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Figure  4.   Theoretical  form  of  a  poikilo thermic  temperature- 
development  profile.   (A)  upper,  nonlinear  temperature 
inhibiting  range.   (B)  middle,  linear  temperature 
dependent  range.   (C)  lower,  nonlinear  temperature 
inhibiting  range.   (D)  critical  point. 
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Above  the  threshold  point,  developmental  speed  increases  with 
increasing  temperature  to  some  maximum.   This  portion  of  the  curve  is 
considered  to  have  an  almost  linear  relationship  (Pearis  1914, 
Andrewartha  and  Birch  1954,  Laudien  1973).   The  slope  of  the  linear 
portion  of  the  temperature  profile  is  of  biological  meaning  by 
depicting  the  relationship  between  rate  of  development  and  temperature. 
Different  species  held  at  the  same  temperature  will  develop  at 
different  rates  (Howe  1962) . 

"Optimum  temperature"  has  a  rather  nebulous  meaning.   Generally, 
this  is  considered  to  be  the  temperature  that  results  in  the  shortest 
development  time.   Optimal  temperatures  may  differ  with  species, 
races,  developmental  state,  or  individual  process  (Laudien  1973). 
Another  definition  of  optimal  temperature  is  that  temperature  at 
which  the  lowest  mortality  is  incurred  (Laudien  1973) .   The  optimum 
temperature  must  be  interpreted  within  the  context  of  systems  being 
described.   Present  investigation  will  focus  on  a  description  of  the 
developmental  process,  per  se  and  will  not  rely  on  previously  used 
definitions  of  "optimum  temperature. "  Of  concern  is  how  the  process 
varies  with  temperature. 

After  the  temperature  causing  maximum  developmental  rate  has  been 
reached,  a  decrease  in  developmental  rate  with  increasing  temperature 
occurs.   Whether  this  slowing  is  due  to  injury  or  enzyme  inhibition 
is  not  completely  clear  (Laudien  19  73).   In  general,  this  decrease 
is  not  steeper  than  the  slope  of  the  increase  to  the  maximum  (Precht 
et  al.  1973). 


19 


Mathematical  Formulae  for  Development 

Janisch  (1932),  Andrewartha  and  Birch  (1954),  and  Precht  et  al. 
(1973)  discuss  several  mathematical  formulae  generally  used  to 
describe  the  temperature-development  profile  for  poikilo therms.       :. 
Each  of  the  authors  discusses  the  theoretical  basis  of  the  formulae,  , 
the  usefulness,  and  the  ease  of  application.   These  several  formulae 
can  be  divided,  roughly,  into  three  groups:   (1)  those  with  theo- 
retical basis  for  shape,  (2)  those  with  empirically  observed  fit, 
and  (3)  the  logistic  equation  (Janisch  1932) .   The  first  group 
includes  the  Arrhenius  equation,  the  reaction  rate-temperature  rule 
(Van't  Hoff  Rule),  and  the  catenary  formula.   The  second  group  contains 
formulae  for  temperature  summation,  the  parabola  used  by  Pearis, 
and  Belehradek's  formula.   Thirdly,  there  is  the  logistic  curve  of 
Verhulst  and  Pearl. 

The  Arrhenius  equation  characterizes  chemical  reactions  by 
"energies  of  activation."  The  natural  logarithm  of  developmental  rate 
will  change  linearly  with  absolute  (  Kelvin)  temperature,  provided 
that  energy  of  activation  requirements  are  met.   Thus,  simple  chemical 
reactions  can  be  characterized  by  respective  energies  of  activation. 
However,  this  method  appears  inappropriate  since  developmental 
processes  are  long  series  of  complex  reactions,  not  simple  reactions. 
Van't  Hoff  Rule  characterized  chemical  reactions  by  the  "Q,f^  rule." 
/Q,Q  implies  a  doubling  in  reaction  rate  for  each  10  C  increase  in 
temperature^/  However,  Q   values  produced  by  this  formula  are  not 
accurate  (Andrewartha  and  Birch  1954) ;  therefore,  the  use  of  this 
formula  usually  is  reserved  for  a  general  look  at  the  linear  portion 
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of  the  curve  (Figure  4) .   The  catenary  formula  of  Janisch  is  the 
combination  of  two  exponential  equations  to  describe  the  acceleration 
and  retardation  of  development.   While  Janisch  (1932)  maintains  that 
the  shape  of  the  curve  has  a  basis  in  theory  (that  the  shape  is  the 
same  as  that  of  various  chemical  reactions),  Laudien  (1973)  classifies 
this  formula  as  empirical.   Regardless  of  classification,  the 
catenary  is  very  difficult  to  utilize. 

The  temperature  summation  as  in  Pearls  (1914)  and  Belehradek 
(Janisch  1932)  claims  no  basis  in  reaction  rate  theory.   Both  formulae 
produce  linear  profiles,  useful  for  predictions,  within  the  middle 
temperature  range.   The  basis  for  the  wide-spread  use  of  this  rule  is 
the  reasonably  good  precision  in  the  linear  region  of  the  temperature 
profile  and,  moreover,  the  extreme  ease  with  which  it  is  applied. 

The  logistic  curve  of  Verhulst  and  Pearl  (Andrewartha  and  Birch 
1954)  has  no  basis  in  reaction  rate  theory.   However,  its  shape 
more  closely  describes  the  temperature  profile  than  does  the  linear 
temperature  summation  method.   This  curve  is  more  difficult  to  apply; 
however,  accuracy  in  description  of  the  temperature-development 
profile  is  increased  (Stinner  et  al.  1974). 

Temperature  Summation  (Day-Degree)  Model 

The  day-degree,  or  temperature  summation  rule,  is  probably  the 
most  widely  used  mathematical  model  for  temperature-development 
relationships  (Sharpe  and  DeMichele  19  77).   Use  of  the  temperature 
summation  rule  in  the  United  States  began  just  after  the  turn  of  the 
century.   Sanderson  (1910)  supplies  a  review  of  the  state-of-the-art 
at  that  time  and  a  general  discussion  of  the  "critical  point"  and 
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"thermal  constant,"  as  well  as  how  these  might  vary  among  species, 
locations,  and  developmental  stages.   Pearis  (1914,  p.  176)  provides 
some  experimental  evidence  of  how  the  summation  could  be  used  and 
states  the  following: 

1)  The  rate  of  velocity  of  insect  development 
is  affected  by  temperature,  and  other  factors 
being  constant,  this  rate  increases  in  direct 
proportion  with  the  increase  in  temperature, 
yytthin  the  normal  limits  of  development. 

2)  The  curve  expressing  the  increase  in  rate  of 
development  is  a  true  hyperbola. 

3)  The  developmental  zero,  or  what  is  styled   •'.' 
the  "critical  point,"  is  at  or  near  the  point 
where  the  reciprocal  curve  for  the  time  factor 
intersects  the  temperature  axis. 

4)  The  thermal  constant  for  an  insect  or  any 
stage  of  an  insect  is  the  constant  for  the 
developmental  curve  for  such  an  insect  or 
stage. 

5)  The  effective  temperature  for  conditions  of 
variable  temperature,  i.e.,  the  ordinary  daily 
variations,  is  higher  than  the  mean  for  the 
period. 

Janisch  (1932)  evaluates  several  mathematical  expressions  for  use 

as  temperature-development  models.   He  notes  that  the  reciprocal  of 

the  parabola  in  general  use  for  day-degree  summation  is  a  straight 

line  and  thus  only  fits  a  portion  of  the  temperature  development 

profile.   Janisch  also  maintains  that  the  hyperbola  has  no  basis  in 

theory  but  simply  is  an  empirically  observed  fit.   Janisch  maintains 

that  an  "S"-shaped  curve,  specifically  the  catenary  curve,  has  a 

theoretical  basis  and  agrees  empirically  with  experimental  evidence. 

Whether  or  not  one  agrees  with  the  views  of  Pearis  or  Janisch,  the 

fact  remains  that  use  of  the  day-degree  model  is  still  very  prevalent. 
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Though  day-degree  models  are  often  constructed  to  explain 
biological  phenomena,  selected  points  on  the  temperature-development 
profile  are  often  deleted.   Also,  these  models  usually  are  not 
validated;  i.e.,  tested  with  real-world  data  not  used  in  model 
construction.   Dunbar  and  Bacon  (1972)  used  a  day-degree  model  to 
compute  developmental  rates  of  three  species  of  Geocoris.   However, 
they  deleted  high  temperatures  (38.3  C)  and  gathered  no  validation 
data.   Hsieh  et  al.  (1974)  constructed  a  day-degree  model  for  the 
alfalfa  weevil,  but  no  validation  data  were  presented.   Thomas  et  al, 
(1974),  working  with  horn  flies,  developed  a  day-degree  model,  again 
deleting  several  upper  and  lower  temperatures.   Chmiel  and  Wilson 
(1979)  constructed  a  day-degree  model  for  the  meadow  spittlebug, 
Philaenus  spumarius  (L.),  in  both  laboratory  and  field.   Although 
there  was  no  significant  difference  between  the  lines  generated  for 
laboratory  and  field  insects,  neither  was  used  for  any  real-world 
prediction. 

Day-degree  models  can  be  used  for  prediction  depending  on  the 
desired  accuracy.   Butler  et  al.  (1974)  maintain  that  in  specific 
situations,  a  simple  day-degree  model  is  adequate  for  analysis  of 
increases  and  decreases  in  abundance  of  natural  insect  populations. 
Hammond  et  al.  (1979)  developed  a  day-degree  model  for  Plathypena 
scabra  (F.)  using  daily  maximum  and  minimum  (max-min)  temperatures  as 
inputs.   Prediction  of  emergence  of  first  generation  adults  in  the 
field  was  calculated  to  be  July  25th.   Actual  emergence  began  on 
July  19th  in  one  field  with  a  large  increase  on  July  24th,   In  another 
field,  first  detection  was  on  July  24th.   However,  in  the  latter  field, 
adults  were  present  in  substantial  numbers  when  first  detected. 
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Sterling  and  Hartstack  (1979)  used  a  day-degree  model  to  predict  the 
spring  emergence  of  cotton  fleahoppers,  Pseudatomoscelis  seriatus 
(Renter) .   A  fit  to  max-min  data  function  was  used  as  input  to 
estimate  hourly  temperatures,  and  the  model  objective  was  predicting 
the  day  of  50%  fleahopper  emergence.   Several  collections  (17)  of 
field  insects  were  made  for  use  as  validation  data.   Four  of  these 
collections  subsequently  were  not  used;  six  were  used  in  construction 
of  the  model;  the  reniaining  7  field  collections  varied  as  follows: 
1,  1,  2,  4,  4,  6,  and  12  days  late. 

Logistic  Equation  Model 

Several  authors  have  utilized  forms  of  the  logistic  equation  for 
more  accurate  estimation  of  developmental  rates  than  the  linear  day- 
degree.   Trips  (1972)  used  the  logistic  equation  to  calculate  rate  of 
development  of  Tororhynchites  brevipalpis  Theobald.   DeVita  (1974) 
used  the  logistic  equation  in  a  deterministic  simulation  model  for 
growth  (increase  in  larval  biomass)  of  Manduca  sexta  (L.).   Validation 
with  historical  data  was  attempted,  and  the  model  fit  data  sets  with 
a  high  degree  of  accuracy  (p  =  0.002).   However,  validation  was  not 
conducted  independently,  as  the  same  data  were  used  to  estimate 
parameters  within  the  model.   An  algorithm  which  utilizes  the  logistic 
equation  was  developed  by  Stinner  et  al.  (1974).   Here,  the  tempera- 
ture profile  is  inverted  after  the  optimal  (fastest  development) 
temperature  occurs.   These  authors  realized  that  the  assumed  symmetry 
of  the  inversion  is  not  correct;  however,  they  maintained  that  because 
the  mortality  rate  in  this  area  of  the  curve  approaches  100%,  the 
error  induced  by  the  assumed  symmetry  is  negligible.   Estimation  of 
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developmental  times  for  the  cabbage  looper,  corn  earworm,  and  the 
spotted  alfalfa  aphid,  Therloaphis  maculata  (Buckton) ,  produced  much 
less  error  than  a  compared  day-degree  model  (35.4%  vs  6%  maximum 
error  for  day-degree  vs  logistic)  (Stinner  et  al.  1974). 

Poikilotherm  Model 

A  stochastic  model  for  simulating  temperature-dependent  develop- 
ment (the  Poikilotherm  Model)  has  been  developed  in  two  phases.   The 
first,  a  reaction  kinetic  model  (Sharpe  and  DeMichele  1977),  is 
based  on  theories  of  enzyme  kinetics  which,  stated  briefly,  are  as 
follows:  Enzyme  systems  of  insect  development  are  very  complex.     These 
systems  can  be  viewed  conceptually  as  a  single  system  with  a  single 
control  enzyme.      Coupled  with  this  concept  are  three  assumptions:   (1) 
the  control  enzyme's  reaction  rate  controls  the  organism's  development, 
(2)  the  development  rate  is  proportional  to  the  concentration  of 
active  control  enzyme  and  its  rate  constants,  and  (3)  the  control 
enzyme  can  exist  in  two  inactive  and  one  active  states.   These 
assumptions  are  depicted  in  Figure  5.   Energy  State  1  is  a 
catalytically  inactive  state  which  predominates  at  low  temperature; 
Energy  State  2  is  an  active  state  which  predominates  at  the  median 
temperature  range;  and  Energy  State  3  is  a  catalytically  inactive 
state  which  predominates  at  high  temperatures.   Transition  among  states 
is  completely  reversible,  except  that  movement  between  States  1  and  3 
is  assumed  to  proceed  through  State  2. 
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From  this  basic  theory,  the  following  equation  for  mean 
development  has  been  derived: 


% 


(0  -  AH  /T)/R 
Te       ^ 


/  1  / 


(AS  -  AH  /T)/R    (AS^  -  AH  /T) 
1  +  e"     "      +e^     ^ 


/R 


where  R^  represents  developmental  rate,  T  represents  temperature 
(  Kelvin),  and  R  represents  the  gas  constant.   Derivation  of  /~1   7 
can  be  found  in  Sharpe  and  DeMichele  (1977).   However,  it  is  important 
to  know  relevance  of  model  parameters.   The  AS  and  AH  represent 
entropy  and  enthalpy  of  transition  between  Energy  States  1  and  2. 
ASj^  and  AH^^  are  the  entropy  and  enthalpy  of  transition  between 
Energy  States  2  and  3.   Note  that  these  constants  are  in  the 
denominator  of  the  equation,  illustrating  the  inverse  relationship 
between  low  and  high  temperatures  and  the  rate  of  development  at 
those  temperatures.   This  is  expected,  since  these  are  the  regions 
within  which  temperature  inhibits  the  developmental  process.   The 
constants  0  and  AH^  represent  the  concentration  of  active  enzyme  and 
enthalpy,  respectively,  of  the  developmental  process  and  are  in  the 
numerator  to  indicate  a  direct  relationship  with  development.   Output 
of  this  model  is  the  estimation  of  6  thermodynamic  constants  which 
describe  the  unique  temperature-development  profile  for  a  given 
organism.   These  constants  then  become  inputs  into  a  stochastic  cohort 
development  model  (Sharpe  et  al.  1977). 

The  cohort  development  model  is  based  on  the  empirical  observation 
that  developmental  times  of  a  cohort  are  distributed  with  a  definite 
skew  toward  the  longer  times.   This  portion  of  the  model  utilizes  the 
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developmental  rate  constants  (equation  J_   1_/)  for  a  particular 
temperature  regime  from  the  reaction  kinetics  model  and  to  predict 
the  distribution  of  development  times  for  some  life  stage  (e.g., 
distribution  of  adult  emergence) .   Details  of  inputs  to  the 
distribution  model  are  described  in  Sharpe  et  al.  (1977).   The 
poikilotherm  model  was  derived  to  predict  (1)  the  beginning  of 
emergence  (of  some  life  stage) ,  (2)  the  shape  of  the  emergence 
distribution,  and  (3)  the  end  of  cohort  emergence.   Because  this  model 
is  relatively  new,  few  examples  of  its  ability  to  perform  are 
available.   Application  of  both  portions  of  the  model  to  historical 
data  can  be  found  in  Sharpe  and  Del-Iichele  (1977)  and  Sharpe  et  al. 
(1977). 

Barfield  et  al.  (1977)  applied  the  poikilotherm  to  the 
hymenopterous  parasite  Bracon  mellitor  Say  where  predictions  of 
developmental  times  for  constant  and  square  wave  temperatures  were 
judged  acceptable.   A  similar  model  was  developed  for  the  fall  army- 
worm,  Spodoptera  frugiperda  (J.E.  Smith)  (Barfield  et  al.  1978). 
Temperature  regimes  tested  included  constant,  square  wave,  "random 
fluctuation"  and  an  average  hourly  temperature  for  a  week  in  July  at 
Gainesville,  Florida.   In  each  case,  the  poikilotherm  model  adequately 
mimicked  observational  data. 

Temperature  Effects  on  Insect  Pathogens 

Studies  concerning  the  effect  of  temperature  upon  insect 
pathogen-host  interactions  have  produced  information  which  indicated 
that  this  relationship  might  fit  the  assumptions  of  the  poikilo thermic 
model.   Moore  (1973),  working  with  fungi,  demonstrated  that 
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temperature  did  not  affect  the  percent  mortality  but  did  affect  the 
time  to  mortality  of  infected  southern  pine  beetles.   A  granulosis 
virus  was  shown  to  exhibit  the  identical  affect  on  the  codling  moth, 
Laspeyresia  pomonella  (L.)  (Sheppard  and  Stairs  1977).   Ignoffo  (1972) 
also  noted  this  while  testing  Bacillus  thuringensis  thuringensis 
Berliner  on  the  pink  bollworm.   However,  this  is  not  always  the  case. 
The  percent  mortality  induced  by  the  Heliothis  NPV  in  corn  earworm 
was  reduced  with  an  increase  in  temperature.   However,  the  experimental 
temperatures  used  were  very  high  (37.7  to  100  C)  (Silverman  and 
Bullock  1968).  Canerday  and  Arant  (1968)  found  that  in  the  NPV  of  the 
cabbage  looper,  higher  doses  reduced  the  time  to  death. 

The  inability  of  an  insect  pathogen  to  infect  its  host  seems 
to  be  prevalent  only  at  higher  temperatures.   Reduced  mortality  of 
the  greater  wax  moth,  Galleria  melloneila  (L.),  subjected  to  high 
temperatures  was  not  due  to  inability  of  the  virus  (Tipula 
irridenscent  virus)  to  penetrate  the  gut  (Tanada  and  Tanabe  1965) . 
These  authors  were  not  able  to  obtain  infection  even  by  injection. 
Thomas  (1959)  tested  the  effect  of  heat  upon  cabbage  looper  and 
corn  earworm  virus  by  drying  the  virus  on  slides  and  subjecting 
these  samples  to  varying  temperatures.   Despite  such  treatment, 
samples  did  not  lose  virulence.   Temperature  can  have  a  therapeutic 
effect  on  virus  infected  organisms.   High  temperatures  could  be  used 
to  "clean"  the  virus  out  of  Lepidoptera  larvae  if  the  larvae  were 
subjected  to  high  temperatures  for  24  to  40  hours  within  24  to  48 
hours  of  infection  (Bullock  1972,  Tanada  and  Chang  1968). 

Pathogen- induced  mortality  in  insects  is  affected  by  temperature 
in  a  variety  of  ways.   Comparison  of  equivalent  insect-pathogen 
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systems  at  different  experimental  temperatures  will  show  differences 
in  the  t-ime   of  mortality  but,  except  for  cases  of  very  extreme 
temperatures,  no  differences  in  the  amount   of  mortality.   Virus 
exposure  to  an  experimental  temperature  before  it  is  administered  to 
a  test  insect  will  produce  mortality,  while  holding  the  insect  at  the 
same  experimental  temperature  during  the  incubation  period  may  not 
result  in  mortality.   Virus  can  be  inhibited  by  temperature  even 
though  artificially  placed  (injection)  directly  at  the  site  of 
replication  (Haemocoel) .   All  of  these  effects  tend  to  indicate  that 
thermal  effects  are  expressed  by  the  replication  of  the  virus,  and 
virus  replication  is  an  enzyme-dependent  process  (Goodheart  1969)  . 

Probit  analysis  has  been  employed  to  explain  relationship  between 
dosage  level  and  mortality,  and  temperature  and  mortality  (Canerday 
and  Arant  1968,  Stairs  1965).   However,  no  attempt  has  been  made  to 
develop  a  model  which  deals  with  the  dynamic  processes  involving 
ingestion,  infection,  replication,  and  mortality  caused  by  insect 
pathogens. 

The  processes  causing  mortality  of  an  insect  by  an  insect  virus 
are  divided  into  two  stages  for  this  discussion.   The  first  stage 
includes  the  events  necessary  to,  but  occurring  before,  actual 
infection.   The  initial  consideration  is  the  concentration  of  virus 
on  the  food  substrate  and  concentration  of  virus  in  the  insect  gut. 
In  the  former  case,  environmental  parameters  such  as  rain,  temperature, 
and  ultraviolet  radiation  will  effect  the  concentration  of  virus 
available  for  ingestion  (Tanada  1971).   The  insect  feeding  rate  and  the 
rate  of  waste  excretion  will  define  the  amount  of  virus  accumulation 
in  the  insect  gut.   Viral  penetration  of  insect  cells  is  then  assumed 
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to  be  a  probability  function  dependent  on  the  concentration  of  virus 
in  the  gut.  The  first  stage  ends  at  the  time  of  infection  of  insect 
cells. 

The  second  stage  of  viral  induced  mortality  deals  with  replication 
of  the  virus  on  a  cellular  level.   Replication  of  the  virus  is  a 
cyclic  process  occurring  in  isolated  cells  or  tissues.   When  an 
individual  cell  is  invaded  by  a  viral  particle,  translation  of  the  host 
nucleic  acid  is  inhibited  at  the  direction  of  viral  nucleic  acid. 
However,  the  host  cellular  machinery  required  for  protein  synthesis  is 
still  functional.   Viral  nucleic  acid  in  essence  outcompetes  the 
host  nucleic  acid  within  the  cell,  and  the  host  cellular  machinery 
begins  to  produce  viral  proteins.   The  eventual  result  of  viral 
replication  is  cell  lysis,  destruction  of  host  cells,  and  liberation 
of  viral  particles  which  infect  other  host  cells  (Smith  1976). 

Synthesis  of  viral  proteins  like  synthesis  of  host  proteins,  will 
proceed  at  rates  largely  determined  by  temperature  as  governed  by  the 
laws  of  enzyme  kinetics.   It  is  this  level  of  viral  infection  that 
the  poikilothermic  model  of  development  (Sharpe  and  DeMichele  1971) 
will  be  used.   The  poikilothermic  model  is  based  on  the  assumption  of 
a  single  developmental  enzyme,  and  that  the  concentration  of  that 
enzyme  defines  the  rate  of  development.   One  might  assume  that  if  viral 
development  proceeds  by  this  same  theory,  the  concentration  of 
developmental  enzjrme  available  to  the  host  would  be  reduced,  thus 
slowing  insect  development.   However,  this  is  not  considered  to  be  the 
case.   Viral  replication  procedes  in  individual  cells,  thus  in  infected 
cells,  the  total  concentration  of  developmental  enzymes  is  converted 
to  viral  development.   Within  uninfected  cells,  insect  development 
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is  assumed  to  proceed  at  rates  unaffected  by  viral  replication  in 
infected  cells. 

The  point  at  which  insect  mortality  occurs  is  very  difficult  to 
quantify  and  probably  varies  with  each  individual.   Infection  is 
generally  considered  to  result  in  eventual  mortality.   Ultimate 
mortality  then  may  be  defined  by  the  previously  mentioned  probability 
function  which  describes  viral  penetration  of  the  insect  gut.   However, 
development  in  uninfected  cells  proceeds  uninhibited,  and  it  can  be 
theorized  that  some  proportion  of  the  insect  population  might  be  able 
to  "outdevelop"  the  virus.   If  this  is  the  case,  it  is  to  be  expected 
that  survivors  will  have  faster  developmental  rates  than  the  average 
uninfected  population. 
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CHAPTER  III 
VELVETS EAN  CATERPILLAR  DEVELOPMENTAL  MODEL 


Methods  and  Materials 

Insect  Rearing 

All  insects  termed  "laboratory  reared"  were  obtained  from  a 
culture  reared  for  24  generations  on  artificial  media  (see  Greene 
et  al.  1973).   Moths  were  held  in  Plexiglass^  cages  (30  cm  X  30  cm  X 
45  cm)  supplied  with  bleached  muslin  cloth  as  an  ovipositional 
substrate  and  a  solution  of  1  g  honey,  0.1  g  ascorbic  acid,  0.1  g 
methyl  para-hydroxybenzoate  and  6  g  sucrose/ 100  ml  water  as  a  food 
source.   VBC  eggs  laid  during  a  single  night  were  collected  the 
following  morning.   Eggs  were  washed  in  a  5%  Clorox  solution, 
neutralized  in  a  10%  sodium  thiosulfate  solution  and  transferred  to 
a  paper  towel.   When  the  eggs  had  dried,  the  toweling  was  cut  into 
sections  and  taped  to  the  inside  of  lids  to  12  oz  Dixie  Fonda^ 
containers  which  contained  artificial  diet.   Containers  were 
transferred  to  the  appropriate  experimental  temperatures. 

Insects  from  naturally  occurring  populations  were  obtained  as 
follows:   Pupae  were  collected  from  the  upper  root  zone  of  soybean 
plants  (ca.  5  cm)  in  a  local  soybean  field,  placed  in  a  Plexiglass 
cage  (30  cm  X  20  cm  X  45  cm),  and  allowed  to  emerge,  mate,  and 
oviposit.   Green  paper  towels  and  soybean  leaflets  were  supplied  as  an 
oviposition  substrate.   A  10%  sucrose  solution  was  provided  as  food 
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for  adults.  The  cage  was  housed  in  a  small  screen  wire  insectory 
located  on  the  edge  of  a  local  soybean  field.  Temperature  information 
was  collected  by  a  recording  hydrothermograph  (Belfort  Instrument 
Company  model  5-594) .  Nightly  egg  lay  was  harvested  the  following 
morning.  Towels  containing  the  VBC  eggs  were  placed  in  1  gal  ice 
cream  cartons  and  allowed  to  hatch.  As  eggs  hatched,  neonate  larvae 
were  placed  into  60  mm  X  15  mm  plastic  disposable  Petri  dishes. 
Soybean  leaflets  were  supplied  as  a  food  source.  Dishes  were 
transferred  to  the  appropriate  experimental  temperatures. 

Constant  Temperature  Experiments 

Experimental  temperatures  of  12.8,  15.6,  18.3,  21.1,  23.9,  26.7, 
29.4,  32.2,  35.0,  37.8,  and  40.6°C  were  used  to  obtain  estimates  of 
developmental  rates.  Laboratory  reared  neonate  larvae  were  placed 
individually  into  1-oz  (30  ml)  plastic  diet  cups  (Thunderbird  Plastics, 
El  Paso,  Texas  model  111)  containing  artificial  medium  and  held  under 
the  appropriate  experimental  temperature  until  adult  eclosion. 
Humidity  was  not  controlled;  however,  relative  humidity  (RH)  was 
always  >^  50%.  A  photoperiod  of  12L:12D  photophaserscotophase  was 
maintained  for  all  temperatures.   Each  larvae  was  observed  once 
every  24  hrs,  and  the  time  from  oviposition  to  1/2"  larva,  1"  larva, 
pupation  and  adult  eclosion  was  recorded.  Time  of  oviposition  was 
assumed  to  be  12  midnight.   The  data  of  any  larval  mortality  also 
was  recorded. 

Validation  Experiments  -  Laboratory 

All  handling  procedures  utilized  in  the  laboratory  validation 
experiments  were  the  same  as  those  of  the  constant  temperature 
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experiments.   Only  the  temperature  regimes  change.   Laboratory  reared 
larvae  were  held  under  a  fluctuating  temperature  regime  which 
simulated  average  hourly  temperatures  for  a  week  in  July  at 
Gainesville,  Florida  (fluctuation  roughly  between  34  and  22°C) . 
Temperatures  were  achieved  through  use  of  a  cam-controlled  environ- 
mental chamber  (Sub-zero  Freezer  Company  model  Lab-400) .   Wild  insects 
were  held  under  one  of  three  experimental  temperature  regimes:   (1) 
Cam-controlled  summer  (Barfield  et  al.  1978),  (2)  21.1°C  constant,  or 
(3)  26.7°C  constant. 

Validation  Experiment  -  Field 

In  1978,  two  sets  of  field  validation  data  were  collected.   The 
first  utilized  laboratory-reared  larvae;  the  second,  larvae  obtained 
from  naturally  occurring  populations.  Except  for  rearing,  both  sets 
of  larvae  were  handled  in  the  same  manner.   Small  larvae  (ca.  2nd 
instar)  were  placed  individually  with  a  small  camelhair  brush  (Red 

Sable  Artist  726-R.U.S.I.  #1)  upon  the  leaves  of  soybean  plants 

R 
contained  in  a  (3.0  mLX1.8mWX  1.8mH)  Lumite  screen  cage 

(Chicopee  Manufacturing  Company  Style  5003800) .   Each  cage  covered 

6.1  row  meters  of  Bragg  variety  soybean  planted  in  90  cm  rows  at  a 

seeding  rate  of  38  seed  per  row  meter.   Cages  were  "cleaned"  by 

visual  checking  and  removal  of  all  insects  from  all  inner  cage 

surfaces,  soybean  plants,  and  ground  to  a  depth  of  5  cm.   Four  cages 

were  used  in  each  test,  and  each  cage  was  supplied  with  150  larvae. 

Temperatures  within  the  cages  were  taken  by  a  recording  hydrothermo- 

graph.   Each  cage  was  checked  daily  for  adults. 
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In  1979,  a  similar  set  of  experiments  utilizing  wild  insects, 
collected  as  larvae  from  Homestead,  Florida,  were  conducted.  However, 
the  presence  of  a  fungal  pathogen,  N^  rileyi,  prevented  completion  of 
the  exp  er  iment  s . 

Results  and  Discussion 

Survivorship 

Survivorship  among  VBC  cohorts  held  at  various  constant  tempera- 
tures changed  very  little  between  21.1  and  32.2°C  (Figure  6). 
Survivorship  decreased  as  temperatures  approached  either  extreme. 
Maximum  survivorship  was  at  23.9  C  and  was  minimum  at  37.8°C.   Egg 
survivorship  was  not  measured  except  that  no  eggs  survived  at 
temperatures  below  12.8  C  or  above  40.6°C.   Greatest  proportion  of 
mortality  was  incurred  during  the  process  of  pupation. 

The  tendency  to  increased  mortality  at  the  pupal  stage  may  be 
a  result  of  increased  physiological  sensitivity  in  the  insect. 
During  the  pupal  stage,  massive  reconstruction  of  tissue  and  growth 
of  muscles  and  wings  take  place.   Signs  of  emminent  mortality  to 
experimental  VBC  larvae  were  shown  by  larval  contraction  and 
abdominal  tanning.   Dark  spots  typical  of  stress  induced  septicemia 
appeared,  followed  closely  by  death  and  putrif ication. 

Mortality  factors  associated  with  the  constant  temperature 
experiment  were  derived  from  two  sources:   (1)  physiological  mortality, 
or  the  inability  of  some  insects  to  cope  with  the  effects  of  tempera- 
ture, and  (2)  handling  mortality.   Handling  mortality  in  these  experi- 
ments was  extremely  small  (<  . 1%) .   Physiological  mortality  might  be 
subdivided  into  two  categories:   (1)  direct  effect  of  temperature  on 
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the  test  animal's  systems  and  (2)  indirect  effect  of  temperature  on 
physiology  which  was  amenable  to  expression  of  facultative  pathogens. 
All  poikilotherms  at  a  given  experimental  temperature  usually  will  not 
complete  the  total  life  cycle.   Furthermore,  mortality  is 
temperature-dependent,  due  to  a  decrease  in  vitality  as  temperature 
departs  from  optimum  (Laudien  1973).   According  to  several  authors, 
(e.g.,  Laudien  1973),  rapidly  developing  individuals  are  long-lived, 
while  those  that  develop  slower  are  short-lived.   Whether  this 
mortality  is  inflected  randomly  or  that  some  individuals  are  "weaker" 
is  not  yet  adequately  defined. 

Temperature- induced  stress  on  experimental  animals  also  may 
effect  processes  which  indirectly  produce  mortality.   Bucher  (1963) 
gives  one  example.   Under  unfavorable  rearing  conditions  (e.g., 
nonoptimal  temperatures) ,  the  gut  wall  becomes  more  permeable  to 
bacteria.   Excessive  numbers  of  bacteria  in  the  gut  are  frequently 
associated  with  septicemia. 

Development 

A  basic  assumption  of  the  poikilothermic  model  is  inhibition  of. 
the  developmental  process  at  temperature  extremes.   Figure  7  depicts 
Arrhenius  curves  for  various  VBC  size  classes.   Arrhenius  plots 
produce  a  linear  range  where  normal  temperature  kinetics  apply  and 
two  (upper  and  lower)  nonlinear  ranges  where  temperature-dependent 
inhibition  of  development  occurs  (Ingraham  1973) .   Figure  7  indicates 
a  reduction  of  developmental  rate  in  the  higher  temperature  range  by 
an  inversion  of  the  curve.   The  extreme  change  in  slope  between  the 
third  and  second  to  lowest  points  should  be  noted.   This  indicates  a 
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disproportionate  reduction  in  developmental  rate  for  a  uniform  change 
in  temperature.   The  VBC  temperature  profile  is  composed  of  three 
regions:   two  temperature- inhibited  portions  (upper  and  lower)  and 
a  linear,  temperature-dependent  portion  (middle) . 

Another  model  assumption  is  that  developmental  rates  should  be 
distributed  symmetrically  about  some  mean.   This  was  shown  to  be  the 
case  with  boll  weevil  and  the  cotton  fleahopper  (Sharpe  et  al.  1977), 
B.  mellitor  (Barf ield  et  al.  1977) ,  and  fall  armyworm  (Barf ield  et  al. 
1978) .   Frequency  distribution  histograms  of  developmental  rates  for 
the  intervals  oviposition  to  1/2"  larva  (OH)  and  oviposition  to  1" 
larva  (01)  are  shown  for  the  upper  nonlinear  temperatures  (Figure  8) , 
the  linear  temperature  range  (Figure  9) ,  and  lower  nonlinear 
temperature  range  (Figure  10) .   These  graphs  depict  the  probability 
distributions  of  developmental  rates  for  several  life  stages. 
Thorough  understanding  of  these  plots  demands  explanation  as  to  why 
developmental  theory  based  on  symmetry  in  rate  distributions  was 
formulated. 

Initial  model  construction  (Sharpe  and  DeMichele  1977)  noted 
that  distributions  of  reported  developmental  times  for  many 
poikilotherms  were  skewed  toward  longer  developmental  times.   The 
inverse  of  this  empirical  distribution  (i.e.,  rate  =  1/time) 
should  produce  the  distribution  for  rates.   This  inverse  distribution 
was  characterized  most  closely  by  the  Quadratic  Probability  Density 
Function  (see  Sharpe  et  al.  1977).   Comparisons  of  the  rate 
distribution  for  VBC  focus  recognition  of  several  experimental  design 
limitations.   Visual  inspection  indicates  that  the  form  of  VBC  rate 
distribution  is  much  closer  to  the  Normal  distribution  than  the 
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velvetbean  caterpillar  reared  on  artificial  diet  and  under 
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26.7,  23.9,  and  21.1°C. 


42 


Uj 
Q: 

ct  o 

^° 
S 


.<0 

0 

CM 

o 

in 

CJ 

.00 

ro 

CM 

% 

CM 

CM 

o 

-••^ 

o 

'J- 

<M 

0 

CM 
CM 

N 

.0 

IT) 

CM 

^^! 

fD 

CO 

^^ 

CO 

lO 

^^ 

— 

3c -^ 

0^ 

--j^ 

■<J- 

.'T 

Uj  -- 

'l- 

K) 

Uj 

-CM 

cvl 

Q 

.0 

.0 
fO 

« 

- 

^ 

oO 

'to 

CO 
CM 

^^ 

o 

-»~ 

O 

CO 
CM 

K) 

CM 

CM 
'10 

<M 

CM 

0 
10 

0 
CJ 

lo            O 

iCt 

in           0 

vo 

t^           m 

pj 

K            in 

CM 

Id 

Q) 

-a 

u 

c 

CO 

CO 

(U 

u 

ro 

u 

00 

CO 

f-H 

iH 

iH 

M 

•H 

0 

fX  14-1 

t-i 

(U 

f^ 

4-) 

3 

to 

0 

0 

x; 

CD 

c 

« 

CO 

(U 

c 

J2 

0 

J-l 

•H 

0) 

4-1 

> 

3 

i-H 

XI 

0) 

•H 

> 

>-i 

4-1 

M-l 

CO 

0 

•H 

Q 

CO 

•W 

V-l 

• 

0 

-a 

X 

0 

0 

•H 

0 

Vj 

(U 

M 

a. 

0 

0 

y-( 

4-1 

0 

cc 

X 

C 

cx 

0 

•H 

Q 

4J 

CM 

3 

.—t 

Xl 

•  • 

•H 

hJ 

Vj 

CN 

4J 

i~H 

W 

•H 

CO 

•0 

^ 

^■^ 

cu 

•i 

TJ 

c 

Q) 

D 

B 

•H 

•a 

4-1 

c 

v-^ 

CO 

(U 

4J 

4J 

OJ 

CO 

•H 

M 

T3 

r-i 

I-H 

CO 

CO 

4-1 

•H 

c 

0 

Q) 

•H 

e 

14-1 

D. 

-H       . 

0 

4-1    U 

T-i 

V40 

<u 

CO  vO 

> 

(U 

C  iJ-i 

p 

0  -1 

3 
60 
•H 


/(i  ISN30      AiniEIVSlOdd 


43 


Quadratic.   Deviation  of  VBC  distributions  from  the  theoretical  form 

may  be  explained  by  recognizing  two  limitations  in  experimental 

design.   The  first  is  non-optimal  temperature  effects.   Excessively 

low  or  high  temperature  may  affect  the  physiological  systems  of  test 

animals  in  quite  different  manners,  many  of  which  are  subtle  and  not 

perceivable  in  controlled  experimentation  (e.g.,  Andrewartha  and 

Birch  1954,  Laudien  1973).   These  effects  were  not  accounted  for  in 

current  experimental  designs;  consequently,  estimates  of  developmental 

rates  at  temperature  extremes  likely  were  modified  by  uncontrollable 

temperature-development  interaction.   Further  complication  arose 

from  the  increase  in  mortality  incurred  as  experimental  temperature 

deviated  from  optimum.   The  second  limitation  was  increased  in 

estimation  of  error  as  samples  were  taken  at  points  temporally 

more  distant  from  time  zero  (oviposition) .   Box  et  al.  (19  78)  discussed 

the  transmission  of  error  through  a  series  of  samples.   For  each 

sample  taken,  moving  away  from  time  zero,  there  is  an  associated  error. 

Each  sample  distribution  is  taken  from  the  previous  distribution  after 

processes  such  as  mortality  have  occurred.   This  results  in  error 

inflation  which  can  be  categorized  by  the  following  models: 

Error  associated  with  oviposition  E  =  E 

to 

Error  associated  with  hatch        E  =  E  +  E^ 

Error  associated  with  0  to  1/2"    E  =  E  +  E,  +  E,  ,^ 

t    o    h    1/2 

Error  associated  with  0  to  1"      E  =  E  +  K  +  E,  ,„  +  E, 

t    o    h    1/2    1 

Error  associated  with  OtoP       E=E+E+E,^+E+E 

t    o    h    1/2    1    P 

Error  associated  with  0  to  A       E  =  E  +  E,  +  E,  ,^  +  E,  +  E„  +  E 

t    o    h     1/2     1    T    A 

where  E^  represents  the  total  error  associated  with  each  event. 
Associated  with  this  error  inflation  is  a  decrease  in  sample  size  due 
to  mortality. 
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The  following  scenario  can  be  used  to  explain  how  these  errors 
are  incurred.   At  all  experimental  temperatures  there  is  an  associated 
mortality,  though  in  some  cases,  it  is  quite  small.   Selectively 
applying  that  mortality  to  the  longer  developers,  simply  because  they 
are  exposed  to  the  environment  longer,  will  have  the  effect  of 
selectively  removing  the  "tail"  from  the  distribution  of  developmental 
times.   When  this  truncated  distribution  is  inverted  to  the  rate 
distribution,  characteristics  of  the  Normal  distribution  will  appear. 
This  is  quite  feasible  at  temperatures  with  high  mortalities.  S  » 

A  further  interaction  of  temperature  and  development  deals  with 
the  relationship  between  growth  stage  duration  and  total  development 
time.   If  all  constant  temperatures  affect  development  in  the  same    *  f^ 
manner,  it  is  to  be  expected  that  a  temperature  induced  change  in 
duration  in  total  development  time  would  result  in  a  proportional 
change  in  each  of  the  growth  stages.   This  hypothesis  can  be  examined 
by  comparing  a  given  size  class  at  each  constant  experimental 
temperature. 

Converting  total  life  span  to  a  physiological  time  scale  (total 
life  span  =  100%)  might  be  expected  to  result  in  a  constant  proportion 

of  total  life  span  being  occupied  by  a  given  life  stage,  regardless  of 

2 
temperature  (P.J.H.  Sharpe  personal  communication).   Table  1  depicts 

the  relative  percent  of  total  development  time  taken  from  oviposition 

to  1/2"  larva  (OH),  1/2"  larva  to  1"  larva  (IH),  1"  larva  to  pupation 


2 
Associate  Professor  Bioengineering,  Biosystems  Research  Division, 

Department  of  Industrial  Engineering,  Texas  A  &  M  University,  College 

Station,  Texas. 
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(IP),  and  pupation  to  adult  eclosion  (PA).   The  proportion  of  total 
development  size  class  appears  constant;  however,  statistical  analysis 
does  not  confirm  this  cursory  appearance.   Variations  associated  with 
proportions  are  not  distributed  normally  and  standard  Analysis  of 
Variance  (ANOVA)  procedures  are  not  adequate  for  testing.   Sokal  and 
Rohlf  (1969)  suggest  the  transformation  by  arc-sine  of  proportion  data 
which  allows  use  of  ANOVA.   Therefore,  calculated  proportion  of  time 
spent  in  any  given  size  class  by  VBC  larvae  was  transformed  by 
arc-sine  and  analyzed  by  linear  regression  with  temperature  as  the 
independent  variable.   Results  of  this  analysis  indicate  that  the 
proportion  of  time  spent  in  the  stages  oviposition  to  1/2"  larva, 
1"  larva  to  pupation,  and  pupation  to  adult  was  not  constant  with 
respect  to  temperature  (a  =  .05);  regression  for  the  stage  1/2"  larva 
to  1"  larva  was  not  significantly  different  from  zero  (a  =  .05). 

Percent  deviation  (standard  deviation/mean  development  rate) 
was  computed  for  each  experimental  temperature.   The  mean  of  these 
deviations  was  taken  for  input  into  the  stochastic  cohort  model. 
Sharpe  et  al.  (1977)  formulated  the  assumption  that  these  deviations 
would  remain  constant  regardless  of  temperature.   Statistical  analysis 
of  this  assumption  was  not  provided.   Table  2  depicts  the  percent  ' 
deviations  for  the  classes  oviposition  to  1/2"  larva  (OH)  and 
oviposition  to  adult  emergence  (OA) .   The  appearance  of  difference 
among  these  present  deviations  with  reference  to  temperature  could 
explain  some  of  the  error  associated  with  the  model's  ability  to 
predict  VBC  developmental  times.   These  values  compare  reasonably  well 
to  percent  deviation  measures  of  boll  weevil  (Sharpe  et  al.  1977)  and 
fall  armyworm  (Barfield  et  al.  1977). 
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Table  2.   Percent  deviation  of  developmental  rates  of  velvetbean 
caterpillar  larva  held  at  various  constant  temperatures 
under  a  12L:12D  photoperiod  and  on  artificial  diet. 


Temperature        Percent  Deviation 

6/R(T) 

OH  OA 


(°C)  6/R(T) 


15.6 
18.3 
21.1 
23.9 
26.7 
29.4 
32.2 
35.0 
37.8 


.08 

.09 

.10 

.08 

.16 

.09 

.15 

.05 

.14 

.08 

.29 

.05 

.16 

.06 

.15 

.05 

.11 

.04 

.14  .06 


Inputs  to  the  first  algorithm  of  the  poikilothermic  model  are 
mean  developmental  rates  and  associated  standard  deviations.   Output 
of  the  model  consists  of  estimates  of  the  six  thermodynamic  constants 
(see  page  26)  which  describe  the  temperature-developmental  profile  for 
the  VBC.   Table  3  displays  developmental  times  and  rates  derived  from 
the  constant  temperature  experiments.   Table  4  lists  the  values  of 
the  thermodynamic  constants  associated  with  VBC  larval  development. 
Estimation  of  values  was  obtained  by  nonlinear  regression  fit  of 
data  gathered  for  VBC  size  classes  (see  Sharpe  and  DeMichele  1977). 
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Table  4. 


Developmental  constants  for  velvetbean  caterpillar  size 
classes  held  at  constant  temperature  on  artificial  diet 
and  12L:12D  photoperiod. 


Size 

Classes 

Constant 

0  to  1/2 

0  to  1" 

0  to  P 

0  to  A 

Units 

0 
HA 
SH 
HH 
SL 
HL 

24.1445 
15305.4 
98.0972 
25963.6 
-47.2525 
-17480.7 

21.8731 
14422.8 
102.762 
27702.6 
-43.6495 
-16169.3 

24.0399 
14028.6 
103.321 
26076.2 
-28.5375 
-13637.4 

29.2604 
18675.6 
75.4515 
18609.2 
-51.8442 
-18749.5 

Time"-"- 

Cal/Mole 

Cal/Mole-°K 

Cal/Mole 

Cal/Mole-°K 

Cal/Mole 

Differences  among  the  developmental  constants  of  each  size  class  , 
are  indicitive  of  the  changes  in  physiological  processes  associated 
with  some  life  stage  change.   The  larger,  more  positive  values  of  0   iv 
and  HA  for  the  0  to  A  class  indicate  slower  or  more  restricted      ~   ./ 
enzymatic  reactions  (see  Morrison  and  Boyd  1974).   This  might  be 
expected  due  to  the  extreme  structural  changes  involved  in  metamorphesis 
to  adult.   Slight  decreases  in  the  values  of  these  constants 

associated  with  0  to  1/2"  larvp  vc  n  .-.,  i"  i 

x/z  xarva  vs  0  to  1   larva  seem  to  be  indicative 

of  only  slight  insect  structural  changes.   For  example,  the 
structural  changes  involved  with  changing  from  a  small  to  a  large 
larva  are  much  simpler  than  changing  from  a  larva  to  a  pupa  or  to  an 
adult.   Size  class  differences  in  the  high  temperature  constants 
(SH  and  HH)  indicate  that  0  to  A  class  is  less  restricted  (i.e., 

small  free  energy)  and  consequently,  more  tolerant  at  high  temperatures. 

The  increase  of  these  values  at  the  0  to  P  stage  may  be  indicative 

of  increased  complexity  of  structural  components  critical  to 

pupation. 


♦  ' 
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Earlier  in  the  discussion,  selective  mortality  was  introduced  as 
one  possible  explanation  of  the  divergence  of  rate  distributions  from 
the  theoretical  form  of  Sharpe  et  al.  (1977).  A  point  of  support  may 
be  found  in  Table  3.   Insect  mortality  evenly  applied  over  the 
distribution  of  developmental  times  for  37.8°C  would  yield  a  mean 
developmental  time  equivalent  to  the  mean  computed  from  a  larger 
sample.   If  the  mortality  were  applied  to  the  "fast  developing" 
insects,  the  mean  would  be  much  lower  than  the  model  would  suggest. 
If,  however,  the  mortality  were  applied  selectively  to  the  "slow 
developing"  insects,  the  mean  would  be  much  higher  and  thus  closer  to 
the  mean  at  35.0°C.   Figure  11  illustrates  this  point.   A  t-test 
comparing  the  mean  developmental  times  from  oviposition  to  adult 
for  35.0  and  37.8°C  yields  no  significant  differences  (a  =  .05). 

This  explanation  does  not  seem  to  hold  for  mortality  at  the  lower 
temperature  extreme.   Visual  inspection  indicates  that  although  large 
mortality  is  incurred  by  insects  held  at  18.3  and  15.6  C,  developmental 
times  do  not  appear  disproportionately  closer  together  than  the  model 
would  predict.   A  simple  explanation  for  this  is  probably  associated 
with  the  breakdown  in  control  of  molting  by  stress  placed  on  the 
organism.   Of  course,  heat  stress  does  not  occur  at  lower  temperatures; 
however,  Okada  (1969)  reports  mortality  at  low  temperatures 
(15.0°C)  due  to  bacterial-induced  septicemia  in  test  of  a  NPV  on  the 
beet  armyworm,  Spodoptera  exigua  (F.).   Also,  limitations  of  experimen- 
tal design  must  be  considered.   Movement  of  insects  from  rearing 
temperature  (26.6  C)  to  experimental  temperatures  of  15.6  and  18.3  C 
was  a  shock  to  the  test  animals.   However,  transfer  from  rearing 
temperature  to  experimental  temperature  involved  only  eggs,  therefore 
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TEMPERATURE 


Figure  11.   Hypothetical  depiction  of  the  placement  of  a  mean 

developmental  rate  acted  upon  by  selected  mortality:   (1) 
mean  with  mortality  applied  to  slower  developers;   (2) 
mean  with  mortality  applied  equally;   (3)  mean  with 
mortality  applied  to  fast  developers. 
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a  large  percentage  of  mortality  was  probably  exhibited  in  the 
embryonic  stage.   Even  so,  surviving  larvae  may  not  have  begun  to  feed 
until  some  adaptation  period  had  elapsed.   Personal  observation 
indicates  that  feeding  by  insects  subjected  to  cool  temperatures  is 
much  slower,  as  indicated  by  production  of  frass,  than  comparably 
aged  insects  at  moderate  temperatures. 

Reduction  of  feeding  would  place  a  drain  on  the  reserved  food 
supply  (fat).   With  reference  to  the  experimental  temperature  18.3  C, 
most  of  the  larval  mortality  was  expressed  in  the  "small"  larvae, 
indicating  that  these  larvae  starved  to  death.   At  other  experimental 
temperatures  (15.6  C) ,  mortality  was  most  prevalent  at  the  pupal  stage, 
which  could  be  attributed  to  the  inability  of  insect  hormonal  systems 
to  function  properly  outside  optimal  temperature  range  (Chapman  1971). 

Data  collected  on  the  VBC  development-temperature  relationship 
conform  to  the  basic  model  assumption.   However,  these  data  were 
collected  on  laboratory  reared  insects  and  a  discrepency  often  exists 
between  the  biological  identities  of  laboratory  reared  vs  naturally 
occurring  insects.   In  order  to  ascertain  differences,  developmental 
comparisons  of  laboratory  and  "wild"  VBC  were  made  under  two  constant 
and  one  variable  temperature  regimes.   Two  cohorts,  one  laboratory 
reared  on  artificial  diet  and  a  second  consisting  of  insects  from 
naturally  occurring  populations,  were  held  at  each  temperature  (Table 
5).   For  each  pair,  a  test  of  equality  of  variance  was  constructed 
(Steel  and  Torrie  1960,  p.  82).   In  all  cases,  variances  were  found 
to  be  unequal.   As  a  result,  t-tests  were  constructed  utilizing  a 
technique  for  the  assumption  of  unequal  variance  (Steel  and  Torrie 
1960,  p.  81).   The  results  obtained  are  mixed.   In  50%  of  the  trials. 
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no  differences  in  developmental  time  were  found,  and  in  cases  where 
differences  were  found,  only  three  contained  differences  of  one  day 
or  more.   In  that  this  model  was  constructed  to  have  a  resolution  of 
no  more  than  one  day  (observations  were  taken  every  24  hours) ,  only 
three  of  the  12  trials  would  have  caused  additional  error.   Thus,  it 
was  concluded  that  the  differences  in  developmental  rates  between 
laboratory  reared  and  "wild"  insects  would  result  in  only  small  errors 
when  employing  this  model  to  simulate  development  of  naturally  .. 
occurring  insects. 

Possible  differences  between  laboratory  and  naturally  occurring 
insects  require  caution  in  model  evaluation.   There  are  two  aspects 
to  this  evaluation.   First,  does  the  model  perform  its  intended 
purpose;   i.e.,  will  this  model  predict  variable  temperature 
development  for  the  VBC?  This  question  can  be  answered  by  laboratory 
experiments.   Second,  is  the  model  generally  adaptable?   If  the 
predictions  are  correct  for  the  laboratory  strain,  will  the  model 
accurately  predict  the  development  of  naturally  occurring  populations 
of  VBC?   It  is  the  latter  that  is  of  concern  to  those  who  must 
utilize  the  model. 

Simulations 

Evaluation  of  the  VBC  developmental  model  will  focus  on  the 
ability  of  the  model  to  fulfill  two  objectives:   (1)  prediction  of 
developmental  times  under  fluctuating  temperature  regimes,  and  (2) 
prediction  of  developmental  times  for  VBC  larvae  in  naturally 
occurring  populations.   This  will  occur  in  7  sequential  stages,  i.e., 
prediction  of  development  times  for  (1)  original   constant  temperature 
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input  data,  (2)  constant  temperature  data  sets  utilizing  wild  insects 
reared  on  soybean  leaflets  collected  independently  of  the  input  data, 
(3)  variable  temperature  development  for  laboratory  insects  on 
artificial  diet  in  growth  chambers,  (4)  variable  temperature 
development  of  insects  from  naturally  occurring  populations  fed 
soybean  leaves  in  growth  chambers,  (5)  laboratory  reared  insects  in 
field  cages  on  soybean,  (6)  insects  from  naturally  occurring 
populations  held  in  field  cages  on  soybean,  and  (7)  development  of 
naturally  occurring  populations  in  soybean  fields. 

Data  on  actual  insect  development  are  discrete,  while 
developmental  model  predictions  are  continuous.   No  statistical  test 
is  available  for  comparison  of  these  two  data  sets;  therefore, 
criteria  for  evaluation  of  the  VBC  model's  predictive  ability  were 
decided  to  be  as  follows:   timing  of  first  emergence,  peak  emergence, 
and  final  emergence  into  a  particular  stage  relative  to  observed 
data.   A  good  "fit"  for  the  timing  of  these  emergences  is  considered 
to  be  +  1  day  of  the  observed  data.  ■  '  ■  - 

Testing  ability  of  any  model  to  predict  the  data  used  for  model 
parameter  estimation  (input)  is  not  validation.   However,  from  a  >  \.;' 
practical  view,  prediction  of  input  data  can  elucidate  gross 
inadequacies  in  the  empirical  data  base.   If  the  model  is  unable  to 
predict  the  data  on  which  the  parameter  estimates  are  based,  one 
cannot  expect  that  the  model  will  generate  accurate  simulations  on 
independent  data  sets.   For  purposes  of  comparison,  predictions  of 
input  are  segregated  into  two  sections  (Figure  12) :   the  size  class 
1/2"  larva  (0  to  1/2")  at  all  experimental  temperatures  and  the 
temperature  26.7°C  for  remaining  size  classes. 
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Figure  12.   Actual  (histogram)  and  simulated  (curve)  times  of 

emergence  for  velvetbean  caterpillar  held  on  artificial 
diet  and  under  a  12L:12D  photoperiod.   (A)  emergence  into 
the  1/2"  larval  size  class  at  various  constant 
temperatures.   (B)  emergence  into  the  1"  larval  size, 
pupal  or  adult  stage  at  26.7  C. 
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Predictions  at  temperatures  in  the  lower,  nonlinear,  region  are 
generally  characterized  by  extremely  wide  (i.e.,  outlying   individuals) 
emergence  intervals  (see  Sharpe  et  al.  1977,  Barfield  et  al.  1978). 
This  effect  is  caused  by  the  very  large  variance  (due  to  high 
mortality  and  small  sample  size)  associated  with  input  means. 
Distribution  of  development  times  for  the  interval  oviposition  to  the 
1/2"  larva  at  15.6°C  illustrates  this  (Figure  12).   Note  that,  aside 
from  the  main  body  of  emergence  times,  there  is  also  a  group  of 
insects  with  much  longer  emergence  times.   Distributions  containing 
outlying  individuals  are  also  found  at  18.3  and  21.1  C. 

Through  the  midrange  temperatures  (21 . 1-32.2°C) ,  simulations 
reflect  the  data  relatively  well.  Generally,  the  0  to  1/2"  simulations 
begin,  peak  and  end  on  the  same  days  as  do  the  input  distributions. 
The  largest  error  in  placement  of  peak  emergence  time  is  at  21.1  and 
26.7°C,  with  the  simulation  peak  appearing  two  days  late  in  both  cases. 
Simulation  terminates  correctly  for  all  cases  in  the  midrange  except 
23.9  and  29.4°C.   However,  these  discrepencies  exclude  only  a  very 
minor  group  of  larvae  which  again  are  outlying  from  the  distribution. 
The  presence  of  these  "outliers"  has  been  noted  with  other  species 
(C.S.  Barfield  personal  communication).   Possibly,  these  organisms 
were  derived  from  progenitors  which  had  reproductive  fitness  enhanced 
by  producing  offspring  which  could  withstand  wide  environmental 
fluctuations  and/or  extremes  in  environment.   Simulations  become  less 
precise  progressing  temporally  through  life  stages.  Largest  error  in 
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simulation  of  onset  of  emergence  was  two  days  found  at  0  to  P. 
The  biggest  error  in  placement  of  the  peak  was  two  days  at  0  to  1",  and 
the  most  inaccurate  simulation  of  final  emergence  time,  an  error  of 
five  days,  was  at  0  to  A.  The  probable  reason  for  the  drop  in 
precision  from  0  to  1/2"  simulations  from  0  to  A  is  the  decrease  in 
sample  size  of  adult  emergence  times  due,  as  described  earlier,  to 
mortality  and  transmission  of  error. 

The  next  step  in  the  sequence  was  to  apply  the  model  to  wild 
insects  again  at  constant  temperatures.   Wild  insects  reared  on  soy- 
bean leaves  were  held  at  21.1°C  (70°F)  (Figure  13)  and  26.7°C  (80°F) 
(Figure  14).   Prediction  of  the  21.1°C  data  for  all  size  classes 
shows  simulated  emergence  into  the  class  well  after  the  empirical 
emergence  has  stopped.   This  also  occurs  at  26.7  C  (Figure  14)  though 
not  to  as  great  an  extent.   Both  of  these  tests  suffered  heavy  loss  , 
in  sample  size  due  to  the  diseases  Nomuraea  rileyi  and  the  VBC-NPV. 
A  total  starting  sample  of  100  individuals  was  reduced  to  56  (44%  loss) 
at  21.1°C.   The  sample  size  at  26.7°C  was  reduced  from  107  individuals 
to  72,  a  34%  loss.   Because  the  effect  of  disease  on  insect  develop- 
ment is  unknown,  any  insect  infected  had  to  be  discarded  even  though 
it  may  have  reached  one  or  more  size  classes.   Comparing  these 
simulations  to  those  of  laboratory  insects  at  these  temperatures 
illustrates  that  the  differences  lie  in  the  "long  developmental  time" 
portion  of  the  curve.   At  21.1  C,  the  0  to  1/2"  distribution  for 
laboratory  larvae  is  on  day  24  and  on  day  18  for  the  wild  larvae,  a 
discrepancy  of  6  days.   The  same  comparison  is  valid  for  all  stages 
at  26.7  C.   With  respect  to  onset  of  emergence,  the  biggest 
difference  is  one  day  (0  to  1/2",  0  to  1",  and  0  to  P)  with  no 
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Figure  13.   Actual  (histogram)  and  simulated  (curve)  emergence  into  the 

1/2"  larval,  1"  larval  size  classes  and  pupal  or  adult  stages 
for  velvetbean  caterpillar  larvae  obtained  from  naturally 
occurring  populations  and  held  at  the  constant  temperature 
21.1  +  1  C  on  soybean  leaflets  and  under  a  12L:12D  photo- 
period. 
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Figure  14.   Actual  (histogram)  and  simulated  (curve)  emergence  into  the 

1/2"  larval,  1"  larval  size  classes  and  pupal  or  adult  stages 
for  velvetbean  caterpillar  larvae  obtained  from  naturally 
occurring  populations  and  held  at  the  constant  temperature 
26.7  +  1  C  on  soybean  leaflets  and  under  a  12L:12D  photoperiod. 
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difference  at  0  to  A.   The  ending  times  are  quite  different,  however, 
vdth  differences  of  2,  4,  7,  and  8  days  for  0  to  1/2",  0  to  1",  0  to  P, 
and  0  to  A,  respectively.   It  might  be  argued  that  these  differences 
are  explained  by  the  influence  of  rearing  on  the  test  insects.  However, 
the  fact  that  almost  no  difference  appears  in  the  timing  of  onset  of 
emergence,  that  the  most  extreme  difference  is  found  at  the  temperature 
(21.1°C)  with  the  highest  mortality  due  to  pathogens  points  to  the 
insect  pathogens  as  the  causative  agent.  Also,  that  the  pathogens  are 
causing  these  differences  by  inflecting  mortality  upon  the  slower 
developing  organisms  is  suggested. 

If  mortality  is  applied  selectively  to  the  distribution  of  host 
emergence  times,  mean  developmental  rate  estimates  will  be  inaccurate. 
Such  errors  could  change  considerably  the  perception  of  how  the  host 
interacts  with  temperature;  i.e.,  actual  developmental  times  could  be 
significantly  different  from  those  measured.   This  interaction  between 
development  and  mortality  presents  a  major  obstacle  to  measurement  of 
development  in  naturally  occurring  populations  where  a  variety  of 
instability  agents  may  be  acting. 

Predictive  ability  of  the  model  variable  temperature  was  tested 
using  both  laboratory  reared  and  wild  insects.   Both  cohorts  were 
tested  at  the  same  time  in  the  same  environmental  chamber.   Figure  15 
depicts  the  actual  (histogram)  and  simulated  (curve)  emergence  time 
into  the  specific  life  stages.   In  general,  onset  of  simulated 
emergence  is  the  same  as  the  actual  with  the  greatest  difference  (two 
days)  at  the  life  stage  0  to  P  for  the  laboratory  reared  insects.   The 
peak  is  properly  placed  (except  for  0  to  1")  for  both  laboratory 
reared  and  the  adult  simulations.   Ending  dates  are  close  (within  +  1 
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Figure  15.   Actual  (histogram)  and  simulated  (curve)  emergence  into 
the  1/2"  larval,  1"  larval  size  classes  and  pupal  or 
adult  stages  for  velvetbean  caterpillar  larvae  taken  from 
naturally  occuring  populations  and  held  on  soybean  leaflets 
or  laboratory  reared  and  held  on  artificial  diet  under  a 
fluctuating  temperature  regime  and  under  a  12L:12D 
photoperiod. 
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day)  for  the  early  stages  (0  to  1/2",  0  to  1",  and  0  to  P)  and  become 
less  accurate  at  the  0  to  A  stage. 

Two  factors  may  have  caused  the  discrepencies  between  simulated 
and  actual  emergence  times.   First,  mortalities  due  to  both  N_^  rileyi 
and  VBC-NPV  were  present  in  both  cohorts.   The  laboratory  cohort  lost 
6%  of  the  sample  size  to  disease  while  the  wild  cohort  lost  A2%.   If 
these  losses  were  applied  to  the  slow  developers,  it  would  cause  an 
unexpected  large  percent  emergence  at  the  early  development  times. 
These  data  were  taken  on  a  daily  basis. 

Results  of  the  first  field  study  in  1978  were  compounded  by 
several  problems.   The  hydrothermograph  used  to  acquire  temperature 
input  was  housed  in  a  shelter  to  protect  it  from  the  frequent  rain 
showers  occurring  in  that  season.   As  a  result,  the  instnunent  was 
kept  in  shade  constantly,  and  this  resulted  in  lower  temperatures 
than  those  to  which  the  insects  were  subjected.   Also,  the 
instrument  failed  to  function  on  several  occasions  resulting  in  a  loss 
of  data.   Simulation  using  these  data  resulted  in  an  "end  of  file" 
error,  indicating  insufficient  information  for  predictions.   It  is 
important  to  note  that  errors  associated  with  malfunction  of  equipment 
and  collection  of  relevant  microclimate  data  constitute  major 
difficulties  in  collecting  adequate  validation  data.   In  an  effort  to 
make  better  use  of  these  data,  two  corrections  were  made:   (1)  missing 
data  were  estimated  by  averaging  the  temperatures  for  three  days  prior 
and  three  days  after  the  missing  data  and  (2)  a  data  set  was  collected 
in  the  same  time  of  year  in  1979  and  used  to  estimate  the  difference 
between  the  hydrothermograph  temperatures  collected  and  those 
temperatures  which  the  insects  experienced.   The  estimation  procedure 
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was  as  follows:   two  thermocouples  were  placed  in  a  field  cage,  one 
in  the  upper  portion  of  the  canopy  which  received  partial  shading, 
and  a  second  in  the  lower  portion  of  the  canopy  which  was  constantly 
shaded.   Hourly  temperatures  were  collected  (Esterline-Angus  PD2064 
microprocessor)  for  five  days  and  the  difference  between  the  two 
temperatures  was  calculated.   This  difference  was  found  to  be  .64°C  dr 
1.2°F.   The  temperature  input  to  the  model  is  in  1°F  increments; 
therefore,  the  hydrothermograph  temperatures  of  1978  were  increased 
1°F.  The  resulting  simulation  (Figure  16)  was  much  closer  (see  page  55) 
to  the  actual  emergence  times,  but  a  discrepancy  of  five  days  was  still 
present.   Whether  this  is  due  to  inaccuracy  in  the  model  or  due  to  the 
inability  to  provide  accurate  temperature  information  is  not  known. 

Results  of  the  second  field  test  of  1978  are  depicted  in  Figure 
17.   The  onset  and  termination  of  adult  eclosion  for  the  simulation 
is  very  close  to  the  observed  days.   The  peak,  however,  is  not  in  the 
proper  place.   Actual  adult  emergence  shows  a  definite  biomodality 
that  the  model  does  not  predict.   Barfield  et  al.  (1978)  also  noted  a 
bimodality  in  fall  armyworm  adult  eclosion  associated  with  widely 
fluctuating  temperatures.   However,  in  that  case,  the  model  accounted 
for  a  portion  of  the  bimodality. 
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CHAPTER  IV 
VELVETBEAN  CATERPILLAR 
NUCLEOPOLYHEDROSIS  VIRUS 
DEVELOPMENTAL  MODEL 


Methods  and  Materials 


Constant  Temperature  Experiments 

Small  discs  of  VBC  artificial  diet  were  cut,  individually  placed 
into  chambers  of  microtiter  multiwell  tissue  plates,  and  surface- 
contaminated  with  virus  dilutions  (3  yl/1  disc),  using  a  micro- 
applicator.   Preparation  of  the  virus  dilutions  was  as  in  Boucias 
et  al.  (1980).   Third  instar  VBC  larvae  were  placed  into  each  well, 
and  the  plates  were  held  at  26.7°C  under  a  photoperiod  of  12L:12D. 
Larvae  were  allowed  to  feed  on  the  contaminated  diet  discs  for  24  hrs. 
Those  larvae  which  had  consumed  the  whole  disc  were  transferred  to 
uncontaminated  diet  (1  oz  plastic  diet  cups,  Thunderbird  Plastics,  El 
Paso,  Texas  model  111)  and  placed  in  the  appropriate  incub.ators. 

Two  series  of  bioassays  were  used  to  determine  the  effect  of 
temperature  on  NPV  in  the  host  VBC.   At  each  experimental  temperature 
(15.6,  18.3,  21.1,  23.9,  26.7,  29. A,  32.2,  35.0,  and  37.8  +1°C),  four 
cohorts  of  larvae  were  fed  one  of  the  following  virus  dosages,  2.49  X 
10^,  1.25  X  10^,  or  2.49  XIO^,  (polyhedral  inclusion  bodies/larva), 
or  no  virus  (control).   These  dosages  represent  ca.  LDqq,  LDcq,  and 
LD2Q  values  for  third  instar  VBC  tested  at  26.7°C  (Boucias  et  al. 
1980).   The  time  to  death  or  pupation  following   inoculation  was 
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recorded.   Constant  temperatures  of  15,  20,  25,  30,  35,  and  40  +  1°C 
were  used  to  obtain  estimates  of  viral  developmental  rate.  At  each 
experimental  temperature,  a  cohort  of  larvae  was  fed  1.06  X  10" 
PIB/larva.   Each  larva  was  checked  daily  and  the  time  from  inoculation 
to  death  was  recorded.   Tissue  smears  prepared  from  the  dead  larvae 
were  examined  using  phase  contrast  optics.  Larvae  showing  definitive 
nucleopolyhedrosis  were  used  in  the  analysis. 

Validation  Experiments 

Two  square  wave  and  two  variable  temperature  regimes  were  used 
to  obtain  validation  information.   In  each  case,  handling  of  test 
animals  was  the  same  as  with  constant  temperature  experiments.   The 
square  wave  temperature  regimes  consisted  of  12:12  hr  at  (A) 
15.6:26.7  +  1°C  or  (B)  21.1:32.2  +  1°C.   The  variable  temperature 
regimes  were  (A)  hourly  temperatures  for  a  week  in  July  at  Gainesville, 
Florida  (fluctuating  roughly  between  35  and  22  C)  and  (B)  a  "random" 
fluctuating  temperature  regime  (fluctuating  roughly  between  15  and 
36°C) .   Temperatures  were  obtained  through  use  of  a  cam-controlled 
environmental  chamber   (Percivil  model  I-35LL) .         -      ■..■■-        .  '  ■* 

Results  and  Discussion 

Development 

The  effect  of  temperature  upon  viral  development,  as  indexed  by 
time  of  death  of  host,  is  of  primary  concern.   Before  this  could  be 
determined,  the  effect  of  dose  level  upon  developmental  time  had  to 
be  measured.   Figure  18  illustrates  the  time  to  death  of  VBC  infected 
with  one  of  three  viral  doses,  at  each  experimental  temperature. 
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Figure  18.   Mean  days  from  inoculation  to  death  for  3rd  instar 

velvetbean  caterpillar  larvae  subjected  to  challenging 
doses  of  velvetbean  caterpillar-nucleopolyhedrosis  virus 
at  various  constant  temperatures  on  artificial  diet  and 
under  a  14L:10D  photoperiod. 
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ANOVA  indicated  that  both  dose  (Table  6)  and  temperature  (Table  7) 
contributed  significantly  (a  =  .05)  to  variation  in  time  from 
inoculation  to  death.   The  time  from  inoculation  to  pupation  for  VBC 
larvae  surviving  one  of  three  challenging  doses  of  VBC-NPV  is 
depicted  in  Figure  19.   Temperature  (Table  8)  and  dose  of  VBC-NPV 
(Table  9)  were  shown  by  ANOVA  to  have  significant  effect  (a  =  .05)  on 
this  time  period. 

The  poikilo thermic  model  was  constructed  to  explain  differences 
in  developmental  times  caused  by  temperature.   The  model,  therefore, 
cannot  account  for  the  differences  in  viral  developmental  times  due 
to  dose  level.   The  previously  documented  experiments  (Tables  8-9) 
illustrate  that  levels  of  viral  dosage  cannot  be  eliminated  as  having 
an  effect  upon  measurement  of  developmental  times  necessary  to 
estimate  parameters  of  the  poikilo thermic  model;  thus,  a  very  large 
(1.06  X  10  PIB/larva)  single  dose  was  used  to  inoculate  test  insects. 
Also,  age  of  test  animals  can  have  an  effect  on  time  from  inoculation 
to  death  (see  Boucias  et  al.  1980);  therefore,  a  single  instar  (3rd) 
was  utilized  in  these  tests. 

In  order  to  discern  how  developmental  systems  of  the  VBC-NPV 
systems  conform  to  the  general  poikilo thermic  theory,  the  same  basic 
assumptions  investigated  for  the  host  (VBC)  must  be  evaluated  for  the 
virus.   The  Arrhenius  curve  for  VBC-NPV  development  from  inoculation 
until  death  is  shown  in  Figure  20.   The  same  form  as  observed  for 
VBC  (Figure  7)  appears  for  the  VBC-NPV.   The  maximum  development  rate 
occurs  at  the  same  temperature  (30-32°C) ,  and  the  inversion  indicating 
temperature  inhibition  occurs  at  temperatures  higher  than  this; 
however,  the  extreme  change  in  shape  noted  at  the  lower  temperatures 


71 


Table  6.  Mean  days  from  inoculation  to  death  for  velvetbean 

caterpillar  larvae  infected  with  nucleopolyhedrosis  virus 
and  held  under  a  14L:10D  photoperiod  on  artificial  diet. 


Dose  N  Mean* 


Control  2  26.5  +  0.7"^  A 

249  PIB/larva  93  14.3+8.8  B 

1240  PIB/larva  127  13.1  +  7.0  C 

2490  PIB/larva  249  11.5  +  6.0  D 


Mean  followed  by  same  letter  not  significantly  different,  DMRT 
a  =  .05. 

+  STD.  DEV. 


Table  7.  Mean  days  from  inoculation  to  death  for  velvetbean 

caterpillar  larvae  infected  with  nucleopolyhedrosis  virus 
and  held  at  various  constant  temperatures  under  a  14L:10D 
photoperiod  on  artificial  diet. 


Temperature  (°C)  N  Mean* 


15.6  76 

18.3  115 

21.1  83 
23.9  113 

26.7  85 

29.4  112 

32.2  102 
35.0  108 


Mean  followed  by  same  letter  not  significantly  different,  DMRT, 
a  =  .05. 

+  STD.  DEV. 


26.8  +  5.5"^ 

A 

22.6  +  4.1 

B 

15.2  +  3.0 

C 

9.9  +  1.4 

D 

9.2  +  1.2 

D 

8.4  +  1.3 

DE 

7.7  +  1.6 

DE 

7.0  +  2.5 

E 
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Table  8.  Mean  days  from  inoculation  until  pupation  for  velvetbean 
caterpillar  larvae  surviving  challenging  doses  of 
velvetbean  caterpillar  polyhedrosis  virus  at  various 
constant  temperatures  under  a  14L: lOD  photoperiod  on 
artificial  diet. 


Temperature  (°C) 


N 


Mean* 


15.6 
18,3 
21.1 
23.9 
26.7 
29.4 
32.2 
35.0 


76 
115 

83 
113 

85 
112 
102 
108 


26.8  +  3.4"^ 

A 

22.6  +  1.5 

B 

15.2  +  1.3 

C 

10.0  +  0.8 

D 

9.2  +  0.7 

E 

8.5  +  1.0 

F 

7.7  +  0.6 

G 

7.0  +  1.0 

H 

*  Mean  followed  by  same  letter  not  significantly  different,  DMRT, 
a  =  .05. 


+  STD.  DEV. 


Table  9.  Mean  days  from  inoculation  until  pupation  for  velvetbean 
caterpillars  surviving  one  of  three  challenging  doses  of 
nucleopolyhedrosis  virus  or  a  control,  reared  on  artificial 
diet  under  a  14L:10D  photoperiod. 


Dose 


N 


Mean* 


Control 
249  PIB/larva 
1240  PIB/larva 
2490  PIB/larva 


311 

232 

186 

65 


13.2  +  7.1"'"  A 
12.9  +  6.5  AB 
12.8  +  7.0  BC 
12.5  +  6.7  C 


Mean  followed  by  same  letter  not  significantly  different,  DMRT, 
a  =  .05. 
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+  STD.  DEV. 
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on  artificial  diet  and  under  a  14L:10D  photoperiod. 
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with  the  VBC  Arrhenius  curve  is  not  evident.   Current  experimental 
design  does  not  allow  investigation  of  this  problem.   Further  work  is 
necessary  to  evaluate  if  this  difference  is  indicative  of  insect 
viruses. 

The  assumption  of  symmetry  with  respect  to  the  distribution  of 
developmental  rates  (Figure  21)  seems  to  be  correct  only  for  those 
distributions  at  midrange  temperatures.   Visual  inspection  indicates 
symmetry  at  25.0  and  30.0  C.   Distributions  at  temperature  extremes 
(15.0,  35.0,  and  37.8  C)  do  not  appear  to  conform  to  this  theory.    ; 

Though  not  a  life  stage,  the  time  from  inoculation  to  death  has 
been  shown  to  be  a  constant,  regardless  of  temperature,  proportion  of  . 
the  life  span  of  survivors  (Canerday  and  Arant  1968) .   The  proportion, 
mean  time  from  inoculation  to  death,  divided  by  time  from  inoculation 
to  pupation  of  survivors  was  computed  for  various  constant  temperatures 
(Table  10) .   Analysis  by  linear  regression  of  arc-sine  transformed  data 
indicates  that  this  proportion  is  not  constant  for  the  VBC-NPV  (a  =  .05) 


Table  10.   Mean  percent  of  a  total  life  span  spent  in  the  viral 

development  stage  of  third  instar  velvetbean  caterpillar 
larvae  infected  with  nucleopolyhedrosis  virus  and  held  at 
various  constant  temperatures  under  a  14L:10D  photoperiod 
on  artificial  diet. 


Temperature 

Days  To 

Days  to 

%  Total  Span 

(°C) 

Death 

Pupation 

In  Virus  Stage 

15.0 

18.1 

^ 

20.0 

8.0 

13.2 

60 

25.0 

6.6 

10.8 

61 

30.0 

5.5 

8.2 

67 

35.0 

5.8 

7.5 

77 

37.2 

6.1 

6.8 

89 

*  No  measurements  taken 
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Figure  21.   Developmental  rate  (time  )  distributions  for  3rd  instar 
velvetbean  caterpillar  larvae  infected  with  velvetbean 
caterpillar-nucleopolyhedrosis  virus  and  held  at  various 
constant  temperatures  on  artificial  diet  and  under  a  14L:10D 
photoperiod. 
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The  percent  deviation  for  VBC-NPV  develonmental  rates  was 
comouted  to  orovide  inout  for  the  stochastic  cohort  model  (Table  ID , 
Like  the  oercent  deviations  for  the  VBC  size  classes  CTable  2) , 
these  do  not  acoear  to  remain  constant  reeardless  of  temnerature. 


Table  11.   Percent  deviation  associated  with  developmental  rates  of 

velvetbean  caterpillar-nucleopolyhedrosis  virus  at  constant 
temperatures  under  a  14L:10D  photoperiod  on  artificial  diet. 

Temperature  (°C)         %  Deviation  /76/R(T)  7 


15.0  .16 

20.0  .12 

25.0  .14 

30.0  .22 

35.0  .22 

37.2  .22 


Although  there  is  some  deviation  from  the  assumptions  of  the  poikilo- 
thermic  model,  data  for  the  VBC-NPV  experiments  seem  to  be  useful  in 
poikilothermic  model  parameter  estimation.   Inputs  into  the  poikilo- 
thermic  model  are  the  mean  rates  and  associated  standard  deviations  for 
the  time  period  "inoculation  to  death"  (Table  12) . 


Table  12.   Developmental  rates  for  velvetbean  caterpillar-nucleopoly- 
hedrosis virus  measured  at  constant  temperatures  under  a 
14L:10D  photoperiod  on  artificial  diet. 


Temperature  (°C)  Rate  (1/T) 


15.0  ,0564  +  .0091* 

20.0  .1263  +  .0155 

25.0  .1534  +  .0225 

30.0  .1901  +  .0434 

35.0  .1705  +  .039 

37.8  .1631  +  .0373 


*  +  STD.  DEV. 
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Output  of  this  model  are  the  six  thermodynamic  constants  which 
describe  the  temperature-developmental  curve  (Table  13).   Like  the 
VBC  model  (Table  4) ,  estimates  of  these  values  were  obtained  by 
fitting  the  equation  for  poikilothermic  development  (p.  24)  to  data 
gathered  at  constant  temperatures  with  a  nonlinear  regression 
technique  (see  Sharpe  and  DeMichele  1977). 


Table  13.   Thermodynamic  constants  obtained  by  fitting  the  equation 
for  poikilothermic  development  to  data  gathered  at 
constant  temperatures  under  a  14L:10D  photoperiod  on 
artificial  diet  for  the  velvetbean  caterpillar- 
nucleopolyhedrosis  virus. 


Constant  Value  Units 


0  33.6773  Time  * 

HA  19027.3  Cal/mole 

SH  94.3536  Cal/mole 

HH  22910.5  Cal/mole 

SL  -45.8154  Cal/mole 

HL  -18137.6  Cal/mole 


*  See  Sharpe  et  al.  (1977) 

The  attempt  to  apply  the  poikilothermic  model  to  the  VBC-NPV 
system  is  complicated  by  the  failure  of  the  data  to  meet  one  of  the  ^ 
model  assumptions.   As  was  noted  with  the  Arrhenius  curve  (Figure  20), 
the  lower  temperatures  do  not  exhibit  the  temperature  inhibition  to  '' 
the  extent  as  was  noted  with  the  hosts.   Consequently,  errors  resulted 
in  estimation  of  the  values  SL  and  HL.   Experiments  are  currently 
on-going  to  obtain  values  of  SL  and  HL  for  temperatures  lower  than 
those  of  this  experiment.   However,  for  the  short  run,  a  fictitious 
point  was  used  which  would  produce  the  correct  theoretical  shape  for 
the  curve,  and  the  constants,  which  appear  in  Table  13,  were  recomputed. 
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The  criticism  that  this  modification  is  simply  an  attempt  to 
"manipulate  numbers"  in  order  to  produce  the  desired  results  is  not 
considered  to  be  the  case  for  the  following  reasons:   (1)  the 
assumption  is  made  that  the  viral  development  does  follow  the 
poikilothermic  theory  and  that  the  lower  values  have  not  yet  been 
accurately  estimated,  and  (2)  the  problem  is  limited  by  inability  of 
the  algorithm  to  fit  the  curve  with  the  current  data.   Obtaining 
accurate  estimates  of  SL  and  HL  for  VBC-NPV  requires  further 
investigation. 

Simulations 

Evaluation  of  simulated  time  to  death,  like  VBC  simulations, 
could  not  be  tested  statistically  because  the  empirical  data  are 
discrete  and  simulation  is  continuous.   Therefore,  the  timing  of  onset, 
peak,  and  final  mortalities  were  compared  visually.   Three  out  fo  four 
simulations  occurred  at  the  square  wave  program  of  15.6/26.7  C 
(Figure  22) .   This  is  to  be  expected  as  the  lower  temperatures  of 
this  regime  reflect  the  inaccuracy  of  the  fictitious  point.   However, 
other  simulations  vary  less  than  two  days  for  any  point.   The  largest 
error  in  starting  time  is  two  days,  and  this  was  associated  with  the 
"random"  (see  Barfield  et  al.  1978)  temperature  regime  (Figure  23). 
Peak  mortality  was  one  day  off  in  the  "random,"  "summer"  and 
21.1/32.2°C  regimes.   Error  associated  with  the  final  mortality  was 
greatest  with  the  "random"  temperature  regime  and  discrepancy  was 
two  days.   In  general,  the  VBC-NPV  model  produced  results  which 
compare  very  closely  with  observed  data. 
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Figure  22.   Actual  (histogram)  and  simulated  (curve)  mortality  of  3rd 

instar  velvetbean  caterpillar-nucleopolyhedrosis  virus  at  a 
square  wave  temperature  regime  of  21.1/23.2°C  or  15.6/26. 7°C 
on  artificial  diet  and  under  a  14L:10D  photoperiod. 
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Figure  23.   Actual  (histogram)  and  simulated  (curve)  mortality  of  3rd 
instar  velvetbean  caterpillar  larvae  infected  with 
velvetbean  caterpillar-nucleopolyhedrosis  virus  at  one 
of  two  fluctuating  temperature  regimes  on  artificial  diet 
and  under  a  14L:10D  photoperiod. 


CHAPTER  V 
CONCLUSIONS 


Studies  were  constructed  to  elucidate  the  effects  of  temperature 
on  development  of  A^^  gemmatalis  and  an  associated  nucleopolyhedrosis 
virus.   Resulting  data  illustrated  a  variety  of  biological  properties 
and  provided  input  into  models  describing  development  of  both 
A.  gemmatalis  and  the  nucleopolyhedrosis  virus.   These  models  will 
be  useful  in  the  formation  of  population  models  of  A_^  gemmatalis . 

Temperature  had  a  significant  effect  on  development  of  both 
A.  gemmatalis  and  the  NPV.   Generally,  developmental  rate  of  both 
systems  increased  with  increasing  temperature  up  to  32.2  C.   At 
higher  temperatures  (35.0  and  37.8  C) ,  developmental  rates  declined. 
Temperatures  between  21.1  and  32.2  C  resulted  in  the  least  mortality 
physiologically  for  A.  gemmatalis.   Extreme  temperatures  15.6,  18.3 
36.0  and  37.8+1  C)  resulted  in  high  percent  mortality,  which  was 
most  pronounced  during  the  metamorphosis  into  the  pupal  stage. 
Virus  replication  was  completely  inhibited  at  40.0  C. 

Present  models  were  based  on  principles  of  poikilo thermic 
development,  as  described  by  general  laws  of  enzyme  kinetics.   General 
success  was  achieved  in  describing  the  observed  data  by  both  models. 
However,  the  developmental  models  were  constructed  to  describe 
responses  to  temperature.   Temperature  is  certainly  a  key  factor 
affecting  development  of  both  systems;  however,  there  are  other  factors 
which  were  not  directly  investigated.   Factors  such  as  age  of  the 
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host  and  dose  level  will  affect  the  developmental  rate  of  the  NPV. 
The  basic  structure  of  thermodynamic  interaction  will  provide  other 
workers  with  a  foundation  upon  which  to  incorporate  the  effects  of 
other  parameters  and  will  allow  the  use  of  sensitivity  analysis  to 
describe  the  resulting  interactions. 

Need  for  the  present  research  was   twofold.   First,  A^   gemmatalis 
is  a  major  pest  of  soybean  in  Florida.   Although  much  information  on 
the  biology  of  this  pest  exists,  insufficient  data  were  available  to 
provide  a  comprehensive  framework  for  estimating  the  temperature- 
dependent  development  of  this  insect.   Elucidation  of  the  basic 
characteristics  governing  the  population  biology  will  allow  meaningful 
evaluation  of  management  tactics.   Disease  induced  mortality 
constitutes  a  heretofore  unquantif ied  variable  in  A_^  gemma talis 
population  regulation.   Description  of  this  process  will  provide 
more  complete  understanding  of  population  biology.   Second,  a  current 
effort  to  construct  an  overall  model  depicting  the  Florida  soybean 
agroecosystem  demands  estimation  of  certain  biological  parameters 
critical  to  the  formulation  of  population  models. 

Results  presented  here  are  a  foundation  but  must  be  recognized 
as  the  first  of  many  steps  necessary  for  understanding  the  VBC-soybean 
system.   Though  some  light  has  been  shed  on  the  development  of 
A.  gemmatalis  and  the  associated  NPV,  much  information  is  still      , " 
lacking.   Areas  of  research  which  need  investigation  concerning 
A.  gemmatalis  development  are  (1)  the  ability  of  A_^  gemmatalis  lairvae 
to  adapt  to  constant  temperatures,  (2)  quantification  of  physiological 
mortality  and  description  of  how  this  mortality  is  applied  to  cohorts 
of  developing  larvae,  and  (3)  effect  of  temperature  on  specific  ages 
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of  larvae.   The  NPV  system  requires  additional  investigation  on  the 
following:   (1)  effect  of  dose  level  on  development  time,  (2)  effect 
of  age  on  development  time,  (3)  effect  of  insect  resistance  on  viral 
development,  and  (4)  description  of  how  NPV  mortality  is  applied 
to  cohorts  of  developing  host  larvae.   The  NPV  developmental  model 
can  provide  adequate  estimation  of  when  host  will  die;  however,  the 
most  logical  step  is  to  quantify  how  many  hosts  will  die  given  a 
specific  set  of  conditions  and  where  in  the  cohort  of  developing 

larvae  will  this  mortality  appear. 

3  ^    . 

Barfield   (personal  communication)  has  provided  a  senario 

depicting  a  series  of  steps  through  which  development  of  these  models 
of  useful  form  should  proceed  (Figure  24) .   The  model  describing 
temperature-dependent  development  for  A^  gemmatalis  has  been  completed 
through  Step  3.   Information  concerning  field  validation  has  been 
obtained;  however,  these  data  are  insufficient  for  comprehensive 
validation.   Certainly,  more  field  trails  are  necessary  to  provide 
confidence  in  the  model  preformance.   Although  further  investigations 
are  necessary  for  evaluation,  the  current  form  of  the  model  can  be 
used  to  provide  estimates  of  developmental  rates  for  naturally 
occurring  populations  of  A^  gemmatalis  (Step  5) .   As  updated  forms  of 
the  model  become  available,  estimates  of  development  time  should  become 
more  accurate. 

Construction  of  the  VBC-NPV  model  has  not  reached  the  level  of 
completeness  of  its  counterpart.   Development  of  the  VBC-NPV  model 
has  proceeded  through  Step  3;  however,  factors  other  than  temperature 
must  be  accounted  for  before  the  VBC-NPV  model  can  be  considered 
complete.   Both  the  age  of  the  host  and  the  dose  level  of  VBC-NPV 
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—  Poikilo thermic  Development  Model 


(Sharpe  and  DeMichele  1977) 


~   Stochastic  Cohort  Model 


(Sharpe  et  al.  1977) 


—  Validation  Laboratory 
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—  Incorporation  into  Pest  Management  Plan 


Figure  24.   Construction  scheme  for  parameter  estimation,  validation, 
and  utilization  of  the  poikilothermic  model  for  insect 
development. 
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received  by  the  host  have  a  direct  effect  upon  the  amount  of  time 
required  for  the  resultant  mortality.   This  will  necessitate  a  series 
of  investigations  designed  in  such  a  manner  as  to  describe  how  these 
parameters  will  affect  the  temperature-dependent  development  of 
VBC-NPV.   One  possible  method  of  conducting  these  experiments  might 
be  to  estimate  the  thermodynamics  constant  for  insects  of  varying 
ages  which  have  been  inoculated  with  a  series  of  doses. 

Realization  of  the  use  of  systems  analysis  stems  from  the  need  to 
decipher  complex  systems.   The  A_^  gemmatalis  and  VBC-NPV  developmental 
models  should  provide  a  base  for  understanding  two  of  the  basic 
parameters  of  insect  population  dynamics  in  the  Florida  soybean 
agroecosystem.   Although  the  current  models  proivde  only  incomplete 
information  on  VBC  development  and  mortality,  the  existence  of  these 
models  will  provide  for  future  workers  to  describe  more  easily  other 
characteristics  which  make  A^  gemmatalis  a  major  soybean  pest  in 
Florida. 
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